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1  3  ABSTRACT  (Maximum  200  words) 

'This  final  report,. ^covering  the  work  done  under  contract  DNA  001-84-C-0015 
during  the  performance  period  of  11  October  1983  to  31  January  1986,  describes  our 
research  conducted  in  close’  coordination  with  the  efforts  at  the  Naval  Research 
Laboratory.  During  this  collaboration,  we>used  and  modified  codes  such  as  MANIC, 
developed  by  Mission  Research  Corporation*-*  to  model  vacuum  diode  regions  of  the 
DNA  superpower  ten  terawatt  generators  such  as  the  DNA/Aurora  14  TW  pulse’  at  the 
Harry  Diamond  Laboratories. 

In  regards  to  the  unpredictable  destruction  of  the  tantalum  anode  .uring  shots, 
it  is  shown  that  when  ions  are  present  in  the  diode  gap,  the  electron  flow  pinches, 

resulting  in  a  high  energy  density  at  the  anode.  This  is  believed  lo  be  a  mechan¬ 

ism  of  the  anode  destruction.  When  ions  are  not  present,  our  models  show  that  the 
flow  does  not  pinch.  In  coupling  of  the  MAGIC  code  to  a  model  > f  a  Plasma  Radia¬ 
tion  Source  it  was  also  found  that  the  electron  flow  outside  >f  the  pinch  is  much 
less  important  than  anticipated  with  the  maximum  current  in  free  electrons  only  a 

few  percent  of  that  through  the  pinch  for  tvpical  plasma  ’adiation  sources  as  those 

used  on  multi-terawatt  generators,  v  •  _ 
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SECTION  1 


INTRODUCTION 


This  final  report,  covering  the  work  done  under  contract  DNA001-84-C-0013 
during  the  performance  period  from  11  OCT  83  to  31  JAN  86,  describes  the  results 
of  research  conducted.  The  first  task  in  the  contract  is  to  develop  an  appropriate 
electromagnetic,  fully  relativistic  particle  in  cell  code  which  realistically  models  the 
vacuum  regions  of  the  DNA  superpower  ten  terawatt  generators.  Sections  2,  3  and 
4  address  this  question. 

One  superpower  generator  is  the  DNA/Aurora  pulser  at  the  Harry  Diamond 
Laboratories.  Section  2  describes  a  diode  which  is  modelled  using  a  modified  version 
of  Mission  Research  Corporation's  MAGIC  code.  It  turns  out  that  this  code,  written 
for  this  type  of  problem  over  the  last  few  years  under  DOE  sponsorship  through 
Sandia  National  Laboratories,  is  able  to  handle  the  Aurora  diode.  The  necessary 
modifications  are  in  the  code  diagnostics,  which  are  insufficiently  flexible  for  our 
purpose.  It  is  shown  that  without  ions  in  the  diode,  the  electron  flow  does  not  pinch 
as  expected,  while  with  ions  the  flow  pinches,  possibly  resulting  in  destruction  of 
the  tantalum  anode. 

A  conclusion  from  this  work  is  that  necessary  PIC  codes  can  evolve  on  the 
basis  of  existing  codes.  Such  codes  include  MAGIC  from  MRC,  C-CUBE  (or  ISIS) 
at  Los  Alamos  National  Laboratories,  MASK  at  SAI,  and  others.  For  all  the  PIC 
code  work  under  this  contract  we  have  used  MAGIC. 

An  integral  part  of  the  statement  of  work  is  to  combine  the  PIC  code  with 
other  codes  that  are  used  by  the  DNA  community,  including  z-pinch  codes.  Section 
3  describes  a  modification  to  MAGIC  that  couples  the  code  to  a  model  of  a  Plasma 
Radiation  Source.  This  research  showed  that  the  electron  flow  outside  the  pinch  is 
much  less  important  than  anticipated.  For  a  plasma  radiation  source  typical  of  the 
ones  used  on  multi-terawatt  generators  the  maximum  current  in  free  electrons  is 
only  a  few  percent  of  the  current  through  the  pinch  (for  the  parameters  selected), 
and  the  maximum  occurs  in  the  early  stage  of  the  pinch,  before  peak  radiation 
production.  This  work  is  described  in  Section  4. 

We  have  also  used  the  MAGIC  code  on  bremsstrahlung  diodes  with  low  induc¬ 
tance.  These  diodes  turn  out  to  be  exceedingly  difficult  to  model.  The  reason  is  the 
widely  disparate  spatial  and  temporal  scales  that  occurs  simultaneously  in  such  a 
diode.  I'nder  these  circumstances  it  is  impossible  to  find  approximations  that,  are 
valid  uniformly  throughout  the  diode  region.  Without  approximations,  however,  a 
full  computation  would  require  excessive  amounts  of  computational  resources.  A 
PIC  computation  on  megavolt  diodes  with  less  than  ~  1  Q  impedance  remains  an 
important  problem  to  solve.  Sustained  research  is  needed  to  resolve  this  problem. 
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Task  2  is  to  continue  the  support  of  the  NRL  Advanced  Concepts  Theory 
program  by  providing  atomic  physics  and  chemistry  models  for  plasma  radiation 
sources.  Sections  5  and  6  describe  the  efforts  done  under  this  task.  Out  of  the 
research  of  leakage  electron  flow  grew  Section  5.  Here,  it  is  suggested  that  some 
alternate  load  configurations,  that  could  be  expected  to  enhance  the  desired  x-ray 
output  from  a  plasma  radiation  source,  be  tried. 

Our  background  in  collisionless  plasma  physics  suggested  that  we  evaluate  the 
need  to  consider  imperfectly  collisional  radiative  plasmas.  Section  6  contains  our 
estimates.  This  work  shows  that  the  hot  spots  in  plasma  radiation  sources,  which  are 
responsible  for  much  of  the  harder  x-rays,  may  be  partly  powered  by  non-maxwellian 
electrons  that  preferentially  excite  higher  atomic  energy  levels  than  are  accessible 
by  a  thermal,  maxwellian  distribution.  In  follow-on  work  it  will  be  necessary  to 
evaluate  this  phenomenon  in  more  detail. 

Task  3  covers  the  transfer  of  working  versions  of  the  codes  developed  under  this 
contract  to  DNA,  and  to  provide  user  guides  to  go  with  the  codes.  All  codes  and 
modifications  to  codes  are  available  at  NRL.  Section  2  contains  a  sample  of  code 
input,  and  Section  3  includes  a  sample  of  a  modification  to  the  MAGIC  code.  The 
comments  throughout  are  sufficent  to  guide  others  to  reproduce  our  work,  and  to 
use  it  in  their  research.  The  code  developed  for  Section  4,  the  Fokker-Planck  code 
FOKKER,  is  available  at  NRL. 


SECTION  2 


AURORA  DIODE  SIMULATIONS 


2.1  INTRODUCTION. 

Every  once  in  a  while  the  anode  in  the  Aurora  flash  x-ray  simulator  sustains 
damage,  necessitating  replacement  of  the  converter  foil  at  a  considerable  expense  in 
downtime  of  the  machine.  The  erratic  occurrence  of  anode  damage  shows  that  '  ere 
is  no  obvious  cause  for  the  damage.  Theoretical,  numerical  and  experimental  studies 
of  anode  damage  on  much  lower  impedance  machines  such  as  Gamble  II  at  the  Naval 
Research  Laboratory  prove  that  ion  emission  followed  by  plasma  emission  from  the 
anode  neutralizes  the  electric  charge  in  the  electron  beam,  after  which  magnetic 
pinching  takes  over.  Damage  occurs  when  the  localized  energy  deposition  exceeds 
what  is  needed  for  melting  or  evaporization. 

The  numerical  simulations  described  in  this  report  suggest  that  also  in  Aurora 
the  anode  damage  is  connected  with  ion  emission  from  the  anode.  Without  ion 
emission,  the  electron  beam  hits  the  anode  in  a  ring  with  approximately  the  cathode 
radius;  when  emission  of  protons  is  included  the  electrons  pinch  down  to  the  center 
of  the  anode,  and  remain  concentrated  there  until  the  end  of  the  pulse.  Pinching 
occurs  only  after  the  ions  have  reached  the  cathode,  and  charge  neutralisation  has 
taken  place  in  the  diode’s  center. 

Pinching  in  the  simulations  is  of  course  not  erratic,  in  contrast  to  the  experi¬ 
ence  at  Aurora.  Hence  the  simulations  do  not  explain  the  anode  damage,  although 
they  were  done  toward  this  purpose;  other  mechanisms  must  be  invoked,  e.g.  vary¬ 
ing  amounts  of  background  gas  and  anode  surface  phenomena.  However,  diode 
modeling  of  Aurora’s  diodes  has  not  previously  been  done,  and  the  results  here  are 
interesting  in  their  own  right  for  what  they  predict  about  the  x-ray  pulse.  Namely, 
the  simulations  suggest  that  the  x-ray  pattern  should  be  ring-shaped  during  the 
first  10  ns.  but  afterwards  the  ring  collapses  on  axis  leading  to  an  intense  central 
spot  in  the  later  part  of  the  pulse. 

Time-integrated  measurements  are  consistent  with  these  results,  however  time- 
resolved  measurements  are  not  yet  available.  The  10  ns  x-ray  risetime  of  the  central 
x-ray  pulse  suggested  by  the  computations  may  motivate  such  measurements.  Dif¬ 
ferences  between  the  measured  data  and  the  simulations  should  give  further  infor¬ 
mation  on  additional  phenomena  left  out  of  the  simulations,  such  as  the  influence 
from  slowly  ionizing  background  gas. 


2.2  SIMULATIONS. 


The  2-1/2  dimensional  electromagnetic  particle  code  MAGIC  used  in  the  simu¬ 
lations  allows  a  non-uniform  grid  shown  in  Fig.  1.  Only  the  front  end  of  the  diode  is 
represented.  The  diode  is  1  meter  long  and  has  a  radius  of  60  cm.  Electron  emission 
is  allowed  from  the  shank  starting  at  20  cm  into  the  diode,  and  from  the  tip  and 
inside  of  the  cathode.  In  the  run  to  be  discussed  ions  are  emitted  from  the  renter 
part  of  the  anode.  Other  runs  had  different  parameters  for  electron  and/or  ion 
emission,  and  various  grid  sizes,  without  appreciable  change  in  the  results.  Section 
2.4  contains  the  input  for  the  run  under  discussion. 
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Figure  1.  Diode  geometry  including  hollow  ring  cathode,  flat 

anode,  and  location  of  the  voltage  and  current  monitors 
at  diode  entrance.  The  grid  is  nonuniform,  decreasing 
from  the  top  inlet  to  the  center  anode  by  a  factor  of  2. 

The  voltage  in  Fig.  2  is  close  to  a  typical  Aurora  voltage,  perhaps  slightly 
shorter.  The  dashed  line  in  Fig.  2  is  the  magnetic  field  Bg  at  the  outer  edge 
(0  60  m)  of  the  diode,  or  the  current  as  determined  from  I{MA)  =  5 Bgr,  where 
Bff(T)  and  r(m)  are  in  MKS  units.  The  voltage  and  current  track  each  other  very 
well,  indicating  a  constant  impedance  of  47  Ohm  in  agreement  with  experience.  It 
may  seem  surprising  that  the  impedance  remains  almost  constant,  even  though  the 
dynamics  of  ions  and  electrons  in  the  gap  changes  considerably  during  the  pulse, 
see  below.  The  reason  is  that  much  of  the  current  is  carried  by  the  electron  sheath 


outside  the  cathode,  where  the  balance  between  electric  and  magnetic  fields  remains 
unchanged  since  no  ions  make  it  to  the  outside  of  the  cathode. 


V 
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Figure  2.  Voltage  (solid)  and  current  (dashed)  at  the  diode  inlet. 

Pinching  of  the  electron  beam  occurs  at  around  40  ns.  Before  the  beam  pinches 
the  electrons  move  as  in  Fig.  3,  at  37.5  ns.  The  electrons  come  off  the  shank 
and  form  a  space  charge  sheath  with  an  edge  that  after  a  small  transition  region 
becomes  roughly  parallel  to  the  cathode.  Electron  emission  from  the  shank  is  not 
quite  uniform,  but  this  is  a  feature  of  the  threshold  field  in  tiie  emission  algorithm. 
Emission  off  the  tip  of  the  cathode  is  uniform,  but  the  electrons  come  oil  in  bunches 
that  propagate  directly  to  the  anode.  The  electron  bunches  reflect  the  emission 
parameters:  electrons  are  emitted  only  after  electric  field  has  had  a  chance  to  build 
up.  From  the  constant  interbunch  distance  in  these  electrons  it  is  evident  that  they 
move  at  the  light  speed.  It  takes  the  electrons  about  3  ns  to  traverse  the  diode  gap. 

The  ions  (protons)  are  shown  in  Fig.  4  at  the  same  time.  At  this  time  the 
ions  have  not  yet  reached  the  cathode.  They  cross  the  diode  gap  rather  uniformly, 
perhaps  partly  attracted  by  the  electron  beam.  The  ions  are  denser  near  the  anode 
because  they  are  emitted  with  a  small  velocity,  hence  they  need  some  time  to 
accelerate  into  the  gap. 

Ion  emission  in  the  simulations  is  not  quite  realistic.  While  field  emission  for 
the  electrons  is  well-established,  including  the  space-charge  limit  coming  from  the 
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Figure  3.  Electron  trajectories  before  pinching,  at  37.5  ns. 

condition  that  the  normal  electric  field  vanish  at  the  conductor  surface,  ion  emission 
is  not  standard.  The  code  uses  the  same  algorithm  for  both  electrons  and  ions,  and 
the  result  is  that  the  ions  are  emitted  almost  simultaneously  with  the  electrons, 
once  the  electric  field  has  exceeded  200  kV/cm.  To  represent  the  suspicion  that 
ions  are  only  emitted  once  the  electrons  have  heated  up  the  anode  surface  the  ion 
emission  is  restricted  here  to  the  center  part  of  the  anode,  somewhat  arbitrarily 
bounded  by  the  impact  radius  of  the  electrons  at  the  time  shown  in  Fig.  3.  ie.  37.5 
ns. 

Electron  and  ion  phase  space  (p:  vs  c)  in  Figs.  5  and  6  corroborate  the  laminar 
development  of  the  electron  and  ion  streams.  The  electron  phase  space.  Fig.  5. 
shows  the  space  charge  sheath  on  the  cathode’s  outside,  between  10  and  40  cm. 
Electrons  there  have  traversed  part  of  the  potential  difference,  about  1  MeV,  and 
are  magnetically  insulated.  Once  they  get  to  the  cathode  edge  the  electrons  pick 
up  energy  as  they  move  toward  the  anode.  Ion  phase  space  is  simple:  each  ion 
behaves  almost  the  same,  simply  accelerating  in  the  applied  potential,  with  a  slight 
variation  in  energy  related  to  the  radial  motion. 

A  short  12.5  ns  later,  at  50  ns,  the  electrons  pinch  as  shown  in  Fig.  7.  The 
space  charge  sheath  on  the  cathode  is  little  modified,  but  as  the  electrons  pass  the 
cathode  tip  they  move  radially  inward  toward  a  point  in  front  of  the  anode  center. 
Their  motion  in  this  region  is  not  properly  represented  in  the  simulation  because  the 
grid  size  in  this  highly  complicated  region  appears  insufficient  to  resolve  the  fields 
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Figure  4.  Ion  trajectories  before  pinching,  at  37.5  ns. 

properly:  the  pinching  coincides  with  the  appearance  of  some  noise  on  voltage  and 
current  (less  than  5  %:  suppressed  in  Fig.  2)  and  the  onset  of  a  loss  in  total  energy 
(up  to  6  %).  Thus  the  details  in  the  pinch  region  are  suspect,  and  will  not  be 
quoted.  However,  the  development  toward  the  pinch  phase  is  properly  computed 
and  should  be  realistic.  There  are  also  a  few  electrons  inside  the  cathode  moving 
away  from  the  anode. 

Ion  trajectories  in  Fig.  8  are  clustered  about  the  electron  beam,  with  continuing 
ion  motion  toward  the  inside  of  the  cathode  and  some  ions  reaching  the  outside  of 
the  cathode  tip. 

As  seen  in  Fig.  0.  the  electron  phase  space  ( p,  vs  c)  still  contains  the  same 
space  charge  sheath  as  12.5  ns  earlier,  but  20  cm  from  the  anode  many  electrons 
lose  their  forward  momentum,  partly  to  be  regained  closer  to  the  anode.  However, 
much  of  the  electron  energy  now  resides  in  radial  electron  motion.  In  addition,  some 
electrons  are  going  backwards  toward  the  inside  of  the  cathode. 

The  ions  inside  the  cathode  still  move  in  a  laminar  fashion,  viz.  the  ion  phase 
space  in  Fig.  10,  but  close  to  the  pinch  the  ions  are  accelerated  beyond  the  diode 
voltage,  apparently  due  to  strong  electric  fields  local  to  the  pinch  region.  Such  de¬ 
tails  will  have  to  be  corroborated  by  better  resolution  in  the  pinch  region.  However, 
it  is  clear  that  the  ions  are  oscillating  about  the  electron  beam,  partly  neutralizing 
its  space  charge  and  allowing  the  pinching  to  occur. 
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Figure  5.  Electron  phase  space:  axial  momentum  pz  versus 
axial  distance  2,  at  37.5  ns. 

Evidently  there  are  many  differences  between  reality  and  the  simulation,  and 
some  of  these  will  be  metioned.  How  ions  are  introduced  into  the  code  appears  par¬ 
ticularly  important.  The  code  assumes  that  ions  are  emitted  whenever  the  electric 
field  exceeds  a  particular  value  (here  2  x  10 7V/m).  While  this  may  be  reasonable 
for  the  field-emitted  electrons  in  reality  the  ions  are  thought  to  be  emitted  when  the 
surface  has  gotten  hot  from  electron  bombardment  followed  by  outgassing,  or  some 
similar  mechanism.  In  any  case  there  is  probably  a  time  delay  in  the  ion  emission 
that  has  not  been  properly  represented  in  the  computation.  A  code  modification  is 
necessary  to  represent  this  feature. 

The  ionic  species  may  also  be  important.  The  computation  assumes  a  single 
species,  singly  charged  hydrogen.  In  reality  there  may  be  a  mixture  of  H+ .  C+ 
and/or  C++ ,  and  probably  other  ions  as  well.  In  any  case  chances  are  that  the 
charge  to  mass  ratio  of  the  actual  mixture  of  ions  is  smaller  that  H+ .  Hence  the  ion 
beam  is  probably  slower  than  assumed  in  the  computation,  and  in  reality  it  should 
take  longer  for  the  pinching  process  to  occur. 

Residual  gas  in  the  diode  gap  may  also  supply  ions.  While  the  electron-ion 
ionisation  cross  sections  are  insufficient  to  account  for  appreciable  ionisation  it  is 
possible  that  ultraviolet  light  generated  by  electron  bombardment  of  the  anode 
ionises  the  background  gas.  Again  a  code  modification  is  needed  to  represent  gradual 
introduction  of  plasma  into  the  diode. 


0) 

s 

E 

o 


Figure  6.  Ion  phase  space:  p ,  vs  z  at  37.5  ns. 

2.3  CONCLUSIONS. 

With  some  insight  gleaned  from  the  simulations  presented  here  it  is  possible 
to  suggest  a  few  remedies  to  the  anode  destruction  problem.  First,  there  is  no 
emission  of  light  ions  if  there  is  no  adhesion  or  absorption  of  these  species  on  the 
anode  surface.  Then  the  electric  field  at  the  anode  will  not  be  able  to  pull  any 
(light)  ions  off  the  surface,  and  there  will  be  no  charge  neutralisation.  While  in 
the  pulse  power  community  a  dirty  vacuum  and  thus  filthy  surfaces  seem  to  be  a 
standard,  a  theoretician  may  rightly  wonder  what  can  be  done  practieall*  with  some 
gettering  in  judicious  places,  etc.  Along  the  same  lines  is  the  suggestion  to  replace 
the  residual  gas  with  a  gas  that  is  hard  to  ionise,  e.g.  helium,  and  to  change  the 
surface  of  the  anode  to  one  that  emits  relatively  little  ultraviolet  light  on  electron 
impact.  Insufficient  information  prevents  selection  of  a  particular  material  here. 

Another  possiblilitv  to  avoid  pinching  could  come  from  a  modification  in  the 
diode  geometry  that  might  keep  the  electrons  away  from  the  axis.  Attempts  at  such 
modifications  have  been  tried  in  the  simulations  without  success. 

It  is  gratifying  that  qualitatively  the  pinch  scenario  is  essentially  the  same  as  in 
low-impedance  diodes.  Without  ions  the  electron  space  charge  prevents  pinching: 
introduction  of  ions  cancels  the  space  charge,  and  allows  magnetic  pinching  to  take 
over.  In  the  present  case,  however,  the  pulse  length  is  comparable  to  the  time  the 
ions  need  to  cross  the  gap.  and  short  compared  to  the  time  any  anode  blowotf  would 
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Figure  7.  Electron  trajectories  after  pinching,  at  50  ns. 

need  to  reach  the  cathode,  at  Acm/fis.  Hence  there  is  no  appreciable  ion  current 
at  the  expense  of  the  electrons,  or  a  quasi-stationary  pinch  state.  Hence  it  appears 
possible  to  avoid  catastrophic  pinching  by  one  or  more  of  the  measures  suggested 
above. 
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Figure  9. 


Figure  10.  Ion  phase  space  pz  vs.  2. 


2.4  CODE  INPUT 


comment  x  aurora  diode  simulation  *  / 

title  x  paur  *  / 

restart  3500  0  / 

timeout  50  1  / 

system  2  / 

fields  1  3 

3500 

0.025e-9 

0.5 

0.5 

0.0 

3 

1.0 

0.2 

0.1111  / 

xlgrid  I  82 

o 

0. 

34 

0.01850 

0.5 

46 

0.01100 

0.5 

x2grid  I  68 

2 

0. 

36 

0.00750 

0.27 

30 

0.00750 

0.30  / 

symmetry  axial 

1 

2,  2 

82,  2  / 

conductor  cathode 

1 

2,38 

36,38  / 

conductor  cathode 

1 

36,38 

36,33  / 

conductor  cathode 

-1 

36,33 

2,33  / 

conductor  cathode 

1 

2,33 

2,2/ 

conductor  anode 

-1 

2,68 

82,68  / 

conductor  anode 

-1 

82,68 

82,2/ 

comment  *  slowly  rising  voltage  * 

/ 

comment  *  mock-up  aurora  pulse 

*/ 

voltage  v(t) 

v(r) 

1.0 

twod 

0.0 

1.0 

1 

2,38 

function  v(r) 

5 

-1 

1.0  / 

function  v(t) 

0 

11 

0.0 

0.0 

10.0e-9 

3.0e6 

20.0e-9 

5.0e6 

30.0e-9 

8.0e6 

40.0e-9 

10.0e6 

50.0e-9 

11.0e6 

60.0e-9 

12.0e6 

70.0e-9 

11.0e6 

80.0e-9 

9.0e6 

90.0e-9 

7.0e6 

100.0e-9 

5.0e6 

particles  emission  null  electrons 

0 

0 

4 

3 

1.75e7  2.5e-7 

1 .25e7 

0  0 

0 

0 

1  8.38 

36.38  / 

particles  emission  null  ions 

1  1 

1 

17 

1  2 

1 

leo  2.5e-5 

2e5 

0  0 

0 

0 

-1  82.42 

82.2  / 

particles  emission  null  electrons 

0  0 

4 

3 

1  2 

1 

1.75e7  2.5e-5 

1.25e7 

0  0 

0 

0 

1  36,38 

36,33  / 

particles  emission  null  electrons 

0  0 

4 

3 

1  2 

1 

1.75e7  2.5e-5 

1.25e7 

0  0 

0 

0 

-1  8,33 

36,33  / 

comment  * 

numerix  data  *  / 

courant 

0  0  / 

kinematics 

1  0 

0 

0 

4/ 

com.  tent  * 

diagnostics  *  / 

trajectory 

250  1 

1.0 

0.0, 1.0 

0.0,0. 6  / 

phasespace 

250  24 

1 

0.0.1.0 

0.0.0. 6 

-5.5e9.5.5e9 

-5.5e9.5.5e9 

-5.5e9.5.5e9  / 

phasespace 

250  31 

1 

0.0.1.0 

0.0.0. 6 

-5.5e9,5.5e9 

-5.5e9.5.5e9 

-5.5e9,5.5e9  / 

phasespace 

250  24 

2 

0.0, 1.0 

0.0.0. 6 

-5.5e7,5.5e7 

-5.5e7,5.5e7 

-5.5e7,5.5e7  / 

phasespace 

250  31 

2 

0.0, 1.0 

0.0,0. 6 

-5.5e7.5.5e7 

-5.5e7.5.5e7 

-5.5e7.5.5e7  / 

comment  * 

voltage  obtained  by 

integrating  el  and  e2  sepa 

ratelv  *  / 

contour 

250  1  2  . 

2 

82.68 

1  2  5 

0 

0  unity  / 

contour 

250  2  2  , 

2 

82.68 

14 


2  > 

5 

0 

0  radius  / 

comment 

*  contour  of  b-theta,  the  current  ’ 

‘/ 

contour 

250 

6  2  . 

2 

82.68 

0  0 

5 

0 

0  actr  / 

comment 

*  contour  of  integral  r  rho  dr  and 

r  jl  dr  *  / 

contour 

250 

10  2  , 

2 

82.68 

2  2 

5 

0 

0  radius  / 

contour 

250 

7  2  , 

2 

82,68 

2  2 

5 

0 

0  radius  / 

function  i 

unity  111. 

/ 

function  radius  5  1  1 

1.0/ 

observe 

10 

0 

twod 

1 

0 

1 

0 

1 

2 

2.38 

2,6S  / 

observe 

10 

0 

twod 

1 

0 

1 

0 

1 

1 

36,38 

82,38  / 

observe 

10 

0 

twod 

1 

0 

1 

0 

1 

6 

68,61 

68,61  / 

observe 

10 

0 

twod 

1 

0 

1 

0 

1 

6 

2,61 

2,61  / 

range 

250 

1 

2 

2,38 

2,68 

0/ 

range 

250 

1 

1 

68,  2 

68,68 

0/ 

range 

250 

1 

2 

68,  2 

68,68 

0/ 

range 

250 

1 

6 

68,  2 

68,68 

0/ 

range 

250 

1 

6 

31,38 

31,68 

0/ 

range 

250 

1 

1 

31,38 

31.68 

0/ 

range 

250 

1 

2 

31.38 

31.68 

0/ 

range 

250 

l 

6 

30,38 

33.68 

0  / 

range 

250 

1 

1 

30.38 

30.68 

0/ 

range 

250 

1 

2 

30.38 

30.68 

0/ 

diagnose 

spacing  0  1 

courant  0  1 

kmatrix  0  1 

pbounds  0  1  / 

statistics 

250  / 

balance 

250 

2.2 

82.68  / 

flux 

250 

1 

36,38 

36.68 

o  / 

flux 

250 

1 

81,2 

81.10 

0  / 

flux 

250 

1 

81.10 

81.20 

0/ 

o 


flux 

25C 

1 

81.20 

81.33 

0  / 

flux 

250 

1 

81,33 

81,40 

0/ 

flux 

250 

-1 

81.2 

81.40 

0  / 

display  1 

2,82 

2.6S  / 

comment  *  begin  the  simulation  *  / 
start  / 
stop  / 
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SECTION  3 


DIODE  CODE  DEVELOPMENT 


3.1  INTRODUCTION. 

Recently  the  electromagnetic  paricle  simulation  code  MAGIC  has  become  avail¬ 
able  to  the  DNA  pulse  power  community.  Over  the  past  several  months  we  have 
become  familiar  with  the  utilization  of  this  code  and  have  developed,  implemented 
and  tested  modifications  that  allow  us  to  begin  modeling  plasma  radiation  source 
(PRS)  diodes.  Our  initial  approach  was  to  model  the  plasma  load  as  an  imposed 
axial  current  density.  This  required  that  MAGIC  be  extended  to  include  in  the  sim¬ 
plest  case  a  stationary  axial  current  density  profile  j:(r)  and  more  appropriately  a 
time- varying  profile  jz(r,t). 

The  first  modification  to  MAGIC  was  to  include  a  stationary,  uniform  axial 
current  density  carried  within  an  annular  ring.  For  this  configuration  the  imposed 
azimuthal  magnetic  field  Bg  is  zero  interior  to  the  ring,  has  a  component  that  is 
directly  proportional  to  r  and  a  component  that  is  inversely  proportional  to  r  exte¬ 
rior  to  the  ring.  Our  approach  was  to  modify  the  Lorentz  force  on  the  particles  to 
include  the  analytic  form  for  this  imposed  Bg.  Section  3.2  presents  our  implementa¬ 
tion  and  tests  of  this  modification.  The  three  cases  described  in  this  section  are:  no 
imposed  current;  +200  kA  current;  and  -200  kA  current.  As  expected  there  is  in¬ 
creased  pinching  in  the  +200  kA  case  and  loss  of  magnetic  insulation  in  the  -200  kA 
case.  The  sudden  appearance  of  a  large  magnetic  field  without  the  corresponding 
inductive  electric  field  is  unrealistic.  Therefore,  implementation  of  a  time-varying 
current  would  be  more  appropriate. 

The  second  modification  to  MAGIC  was  to  include  a  time- varying  axial  current 
density  profile.  Our  approach  in  this  modification  was  to  add  the  imposed  spacially 
and  temporally  varying  current  density  to  that  calculated  from  the  particles.  Con¬ 
sequently.  proper  inclusion  of  the  induced  eiectric  field  would  be  handled  by  the 
existing  field  solvers  in  MAGIC.  Section  3.3  presents  our  implementation  and  test 
of  this  modification.  For  this  case,  we  chose  a  current  that  rises  to  +200  kA  during 
the  first  5  ns  and  remains  constant  thereafter.  As  expected,  the  particle  trajectories 
differ  from  the  constant  current  case,  but  are  qualitatively  similar. 

The  third  modification  to  MAGIC  was  to  allow  the  deletion  of  diagnotics  in 
a  RESTART.  Previously,  it  was  only  possible  to  add  to  the  diagnostics,  rids  cre¬ 
ated  a  situation  where  diagnostics  that  were  no  longer  ol  value  were  required  to 
continue  through  all  subsequent  RESTARTS.  Our  modification  was  designt. 
maintain  the  existing  structure  of  MAGIC  dataset  handling.  Section  3.4  presents 
our  implementation  and  an  example  utilizing  this  modification. 


17 


3.2  MODIFICATION  I. 


A.  Constant  Axial  Current 

A  stationary  current  I  along  the  axis  of  a  cylindrical  diode  as  in  Fig.  1 1  produces 
an  additional  magnetic  field.  The  code  can  include  the  field  as  the  initial  value  lor 
the  magnetic  field  to  be  read  in,  but  it  is  simpler  to  modify  only  the  Lorentz  force 
on  the  particles.  Moreover,  the  analytic  form  of  the  magnetic  field  allows  accurate 
magnetic  gradients  for  future  use  in  drift  orbit  particle  pushing. 


Z 


r 


a.  z-r  view  of  axial  current  ring 


b.  Cross  section  of  current  ring 


Figure  11.  Geometry  for  axial  current. 


The  axial  Z-pinch  is  modeled  by  a  ring  with  constant  current  density 

if  r,  <  r  <  r 

J "  l  0.  elsewhere. 


(1 


Ampere's  law  is 


Figure  12.  Initial  state  without  axial  current  Ic. 
with  bondarv  condition  Bg  =  0  at  the  center  r  =  0.  Then 

Bg  =  0  for  r  <  r;,  (3a) 

Bg  =  ~(>  -  rf/r)/(rj  -  rf)  for  r,  <  r  <  .  (3Ai 

Z7T 

and 

Be  =  7^-  for  r  >  r  j .  (3c) 

l~r 

The  magnetic  field  outside  the  diode  depends  on  the  diode  boundary,  open 
return  current  structures  or  a  closed  cylindrical  shell,  but  no  interesting  particles 
penetrate  to  these  regions  and  this  question  will  not  be  addressed  here,  except  to 
point  out  that  the  typical  ^-dependent  return  current  configuration  is  felt  inside  the 
diode. 

Equations  (3)  are  implemented  in  the  PIC  code  MAGIC  by  a  subroutine 
CCREXT  that  is  called  before  the  particles  are  moved.  Code  modifications  are 
shown  in  Section  3.2. B. 


1!) 


Figures  12,  13  and  14  show  computational  output  solely  as  an  example.  An 
Aurora  diode.  80  cm  long,  50  cm  outer  radius,  connected  to  a  120  cm  long  trans¬ 
mission  line  to  a  voltage  source  rising  in  10  ns  from  0  to  8  MeV.  gives  electron 
trajectories  in  Fig.  12  as  expected.  At  time  t  =  20  ns.  the  particles  suddenly  see 
the  additional  magnetic  field  from  Eq.  (3)  with  current  200  and  -200  kA;  the  result 
is  expected,  namely  increased  pinching  in  the  200  kA  case,  and  loss  of  magnetic  in¬ 
sulation  in  the  -200  kA  case.  Of  course,  the  sudden  appearance  of  a  large  magnetic 
field  without  the  corresponding  inductive  electric  field  is  impossible;  figures  12.  13 
and  14  are  intended  to  suggest  that  axial  external  currents  can  be  included  in  the 
PIC  code,  and  will  produce  interesting  phenomena. 


Figure  13.  Ic  =  200&.4. 
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B.  Listing  of  code  modification 


subroutine  curext 

c  ***********  axial  current  gives  magnetic  field  ********** 
c 

c  provides  a  magnetic  field  due  to  a  constant  axial  current 
c  density  in  a  ring  with  inner  radius  rinc  and  outer  radius 
c  rfic  with  total  current  cure,  the  magnetic  field  is  in  the 
c  theta  direction;  hence,  the  routine  works  in  the  z,r, theta 
c  system  (isym=2)  only,  where  b-theta  =  b3. 
c 

c  the  constant  external  field  bext  should  logically  be 
c  implemented  in  the  present  routine  rather  than  in  the 
c  routine  kinemat. 
c 

c  ***********  variable  definition  ************************* 
c 

c  cure  -  axial  current  (amperes) 
c  rinc  -  mn*»r  radius  of  current  ring  (meters) 
c  rfic  -  find  radius  of  current  ring  (meters) 
c 

*ca, common 
cdirS  nolist 
c 

c  only  perform  for  cylindrical  (z,r, theta)  system 
c  when  cure  is  non-zero, 
c 

c  duplicate  the  next  coding  in  subroutine  contour  for  plotting, 
if  (isvs  .ne.  2)  goto  190 
if  (cure  .eq.  0.0)  goto  190 
c 

cmuO  =  cure*xmuO/topi 
abri  =  abs(rinc) 
rfis  =  rfic**2  -  rinc**2 
c 

do  110  np=l,npgp 
rcur  =  x2pi(np) 
if  (rcur  .le.  abri)  goto  110 
if  (rcur  .ge.  rfic)  goto  120 

b.‘5wf(np)  =  b3wf(np)+cmu0)'frcur-rinc“,‘2/rcur)/rfis 
goto  1 10 

120  b3wf(np)  =  b3wf(np)+cmu0/rcur 

110  continue 


190 


continue 


c  add  constant  external  b-field  here  (to  be  implemented) 
c 

return 

end 


2:5 


3.3  MODIFICATION  II. 


.3.  Time  —  Varying  Axial  Current 


In  contrast  to  a  stationary  current  accompanied  by  a  stationary  magnetic  field, 
a  time-varying  current  gives  an  inductive  electric  field  according  to 


dE 

dt 


=  c2V  x  B  -  — : 

Co 


the  magnetic  field  develops  from 


(4a) 


dB 

dt 


-V  x  E. 


(46) 


and  for  a  current  that  becomes  stationary  after  some  time  the  final  magnetic  field 
should  be  stationary  and  balance  the  j-term  in  Eq.  (4a). 


The  current  density  j  in  Eq.  (4a)  in  the  unmodified  code  comes  from  the 
available  particles.  The  modification  is  to  add.  in  subroutine  JEXTERN.  an  axial 
current  density 


jz{r,t)  =  j{r)I{t)  (5) 

to  the  particle  current  density;  j{r)  is  the  radial  current  profile,  determined  by  a 
function  statement,  and  I(t)  is  the  time  variation.  The  radial  profile  is  normalized 
such  that 


J  j(r)dvT  =  1.  (6) 

In  Eq.  (6)  dvT  is  the  surface  element,  d(wr2),  in  the  radial  direction  (X2  for  isvs 
=  2)  that  is  available  in  the  code  (l/vli2).  The  function  I{t)  can  then  be  given  in 
amperes. 

The  actual  coding  is  attached.  Not  that  in  the  code  the  variable  ji  is  propor¬ 
tional  to  the  current  density  but  includes  various  factors.  In  subroutine  .JEXTERN. 
jt  =  current  density  x  volume  x  timestep  =  charge  displacement  in  a  timestep.  In 
subroutine  JBOUNDS,  j;  is  turned  into  ji  =  -  current  density  x  ( dt/eQ ),  in  addition 
to  enforcement  of  the  boundary  conditions. 

Figures  15  and  16  give  the  output  for  a  sample  problem.  In  this  case  the  current 
rises  to  200  kA  during  the  first  5  ns.  and  remains  constant  thereafter.  The  particle 
trajectories  (Fig.  15)  differ  from  the  constant  current  case  but  are  qualitatively 
similar,  in  particular  the  large  vortex  between  anode  and  cathode  tip.  The  magnetic 
field  (Be)  contours  (Fig.  16)  are  as  they  should  be.  bunched  about  the  axis,  and 
straight  except,  close  to  the  particles.  Contours  2  and  3  give  the  correct  magnetic 
field  (assumed  stationary)  of  0.44T  and  0.7T  respectively. 
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Figure  15.  Electron  trajectories.  Magic  simulation  at  time  1.50E-08  seconds 
for  4  time  steps. 

Potential  inaccuracies  in  the  simulation  are  suggested  by  the  many  wiggles  in 
the  contours,  and  by  the  convergence  of  the  contours  to  the  axis  at  the  right  edge. 
This  last  feature  is  due  to  the  perfectly  conductive  boundary  at  that  edge,  which 
wants  to  put  an  image  current,  on  the  other  side  that  cancels  the  applied  cut  tent. 
Other  boundary  conditions  exist  in  the  code  that  are  more  appropriate. 

The  test  problem  is  sufficiently  unphysical  that  we  should  not  worry  aboui  the 
physical  details,  such  as  boundary  conditions,  wiggles,  and  such.  We  conclude  that 
the  additional  axial  current  is  included  properly. 


Figure  16.  Be  contours.  Magic  simulation  at  time  1.50E-08  seconds 
for  contour  intervals  of  2.70E-01  (MI\S). 


B.  Listing  of  code  modi fication 


subroutine  jextern 


*  external  current  density  j(r.t)  along  axis 
gives  an  external  current  density  profile  (dependent  on  r=x2) 
changing  in  time  along  the  z-axis  in  a  cylindrical 

z,r, theta  (isys  =  2) 

system  only,  radial  and  time  dependence  in  through 
aj  lx  a  j  1 1 

the  current  density  is  normalized  to  this  total 
current  at  each  time  step. 


*****  variable  definitions 


curf  -  total  current  in  current  density  profile/normalization 
jl  -  current  density  *  cell  volume  *  dt 

curt  -  applied  axial  current  dependence  on  time  (amperes) 
cdxr  -  radial  dependence  of  the  current  density 


*ca.common 
cdirS  nolist 

c 

c 

c  the  first  time  around,  initialize  the  current  density 
c  array  cdxr.  and  normalize  to  1  ampere 
if  (isys  ,ne.  2)  goto  'JUO 
if  (nctx  .eq.  0)  goto  900 
if  (curf  .ne.  0.)  goto  200 


c 

do  120  i2=l.i2mx 
call  functio(ajlx.x2hg(i2),cdxr(i2)) 
120  continue 


c  find  total  current  represented  by  given  profile 
do  190  i2= 1  ,i2mx 
curf  =  curf-*-cdxr(i2)/vli2(i2) 
continue 

if  (curf  .ne.  0.)  curf  =  l./curf 
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c 


c  set  external  current  density  in  radius  and  time 
200  continue 

call  functio(ajlt,trun,curt) 
temp  =  -curf*curt*dtimxftaj/epsO 
c 

1  =  1 

do  210  i2=l,i2mx 
do  210  il  =  l,ilmx 

j  1  ( 1)  =  jl(l)+temp*cdxr(i2)/(vlil(il)*vli2(i2)) 

1  =  1+1 

210  continue 

c 

c  jl  contains  current  density  *  cell  volume  *  timestep 
c  jl  is  turned  into  a  source  term  for  curl  e  in  current 
c  boundary  routine,  next  step, 
c 

900  continue 

c 

return 

end 


3.4  MODIFICATION  III. 


A.  Modified  Diagnostics  in  a  RESTART 


In  the  present  version  of  MAGIC  it  is  not  possible  to  modify  all  input  pa¬ 
rameters  and  diagnostics  in  a  RESTART.  However,  it  is  possible  to  add  to  the 
diagnostics,  and  the  result  is  much  more  output  than  desired.  To  avoid  this  prob¬ 
lem,  we  have  added  a  DELETE  feature  that  allows  the  code  user  to  delete  certain 
input  functions  and  output  diagnostics.  MAGIC  counts  the  number  of  these  func¬ 
tions  and  diagnostics,  in  subroutine  DATASET,  and  DELETE  simply  reduces  the 
count  by  the  desired  number. 

The  dataset  definition  is 

DELETE  ADEL  NDEL 

ADEL  -  alphanumeric  label  of  dataset 

to  be  deleted. 

NDEL  -  number  of  these  deleted  datasets. 

The  coding,  a  block  in  subroutine  DATASET,  is  attached.  For  example,  in  the 
original  run  you  have  specified  3  contour  plots;  on  RESTART,  you  want  to  change 
the  contour  interval  on  the  second  plot,  and  the  third  plot  is  useless  so  you  want  to 
delete  it.  To  the  input  data  on  RESTART  you  add 

DELETE  CONTOUR  2 

and  after  this  element  you  add 

CONTOUR  (modified  parameters). 

The  diagnostics  on  RESTART  will  contain  the  first  contour,  and  the  contour  plot 
specified  in  the  RESTART;  the  third  will  have  disappeared. 


Elements  that  can  be  deleted  by  the  present  MAGIC  modification  are: 


COMPLOT 

CONTOUR 

DIAGNOSE 


LINPRINT 

LOGPRINT 

OBSERVER 


DISPLAY  PIIASESPACE 

FLUX  RANGE 

FUNCTION 


As  an  example,  in  Fig.  17  are  the  electron  trajectories  after  a  restart  with  a 
changed  time  dependence  of  the  axial  current,  going  from  200  kA  to  -200  kA  in  5 
ns  between  20  and  25  ns.  The  result  is  intermediate  between  the  200  kA  case  ( 1 ,2b ) 
and  the  -200  kA  case  (1. 2c),  as  expected. 
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B.  Listing  of  code  modification 


c 

subroutine  dataset 


c 

c  ************  element  delete  *********** 
c 

c  (data  sequence) 

c 

c  delete  adel  ndel 

c 

c  (variable  definitions) 

c 

c  delete  -  element  label  for  deleting  datasets  on  a  restart 
c  adel  -  alphanumeric  label  specifying  what  to  delete 
c  delchap  character  variable  containing  adel 
c  ndel  -  number  of  'adel’s  deleted 
c 

c  as  an  example,  if  one  diagnostic  turns  out  to  be  useless 
c  for  your  run  it  can  be  deleted  on  restart;  if  the 
c  diagnostic  is  the  middle  of  three  in  the  initial 
c  specification,  the  last  two  must  be  deleted,  and  the  last 
c  one  must  be  resupplied,  to  modify  a  diagnostic  you  must 
c  delete  it,  and  then  resupply  a  new  diagnostic  in  its  place, 
c 


else  if  (valchar  .eq.  'delete')  then 
nv  =  i 

wi'te(delchar,997)  valu(  1+nv) 
ndel  =  valu(2+nv) 
ndel  =  maxO(O.ndel) 
nv  =  nv+2 

if  (delchar.eq.'complot")  then 

ncpl  =  maxO(O.ncpl-ndel) 
goto4!)l 

else  if  (delchar  .eq.  'contour')  then 
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nctr  =  maxO(0,nctr-nciei) 
goto491 

else  if  (delohar  .eq.  'diagnose')  then 
ados  —  maxO(O.ndos-ndel) 
goto491 

else  if  (delchar  .eq.  'display')  then 
ndis  =  maxO(O.ndis-ndel) 
goto491 

else  if  (delchar  .eq.  'dux’)  then 

nflx  =  max0(0,nfix-ndel) 
goto491 

else  if  (delchar  .eq.  ’function")  then 
nfun  =  max0(0,nfun-ndei) 
goto491 

else  if  (delchar  .eq.  "Imprint")  then 
nlin  =  maxO(O.nlin-ndel) 
goto491 

else  if  (delchar  .eq.  dogprint’)  then 
nlog  =  max0(0,nlog-ndel) 
goto491 

else  if  (delchar  .eq.  ’observe’)  then 
nobs  =  max0(0,nobs-ndel) 
goto491 

else  if  (delchar  .eq.  ’phasespace')  then 
npha  =  maxO(O.npha-ndel) 
goto491 

else  if  (delchar  .eq.  ’range’)  then 

nrng  =  max0(0,nrng-ndei) 
goto491 

c  trying  to  delete  invalid  function 
else 

c 

if  (lsup(isup)  .ne.  1 )  then 
write(6.120)  ndel. delchar 
goto  200 
end  if 
endif 

491  continue 

c 

c  function  deleted:  write  name 
ifflsuplisupl.ne.l )  then 
write(6.1‘21 )  ndel. delchar 
endif 
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call  datasum(0.nv,nval, 0,0,0) 


SECTION  4 


LEAKAGE  CURRENTS  OUTSIDE  AN  IMPLODING  Z-PINCH 

4.1  INTRODUCTION. 

It  is  important  to  know  where  the  current  flows  in  a  radiating  z-pinch1.  Ex¬ 
perimentally,  the  emission  of  kilovolt  x-rays  is  roughly  proportional  to  the  fourth 
power2,3  of  the  current,  whence  a  small  shunt  current  could  have  a  surprisingly 
large  and  deleterious  effect  on  the  pinch’s  radiation  output.  A  quantitative  esti¬ 
mate  of  the  shunt  current  is  the  purpose  of  this  paper.  The  current  outside  the 
z-pinch  plasma  is  generally  thought  to  be  negligible.  For  relevant  parameters  (e.g. 
I  ~  1  M  A  and  pinch  radius  r  ~  1  cm)  there  is  a  large  magnetic  field  in  the  vacuum 
outside  the  pinch  (B  ~  10  T),  and  this  5#  field  should  provide  magnetic  insulation. 
As  shown  below,  however,  magnetic  insulation  in  a  cylindrical  z-pinch  diode  is  not 
necessarily  complete  under  all  circumstances.  Whether  leakage  currents  affect  the 
pinch’s  radiation  output  depends  primarily  on  the  amount  of  material  relative  to 
that  in  the  pinch,  that  is,  the  initial  conditions  in  a  particular  experiment. 

In  the  initial  phase  of  the  electrical  pulse  the  current  through  the  pinch  is 
small  while  the  voltage  is  already  appreciable.  Under  these  circumstances  there  is 
a  transient  electron  current  analogous  to  the  uninsulated  front  of  a  magnetic  insu¬ 
lation  wave.  Thereafter  the  current  through  the  inductive  load  becomes  sufficiently 
large  to  cut  off  this  transient.  Computer  simulation  indicates  that  the  magnetically 
captured  electron  sheath  rolls  up  into  relatively  stable  vortices,  see  Section  4.2. 

After  this  initial  phase  but  still  in  the  early  stage  of  the  pulse  the  z-pinch 
material  is  rapidly  ionized,  a  complicated  process  outside  the  scope  of  this  work.  The 
current  compresses  and  implodes  the  bulk  of  the  plasma,  but  almost  by  definition  the 
current  density  is  low  outside  the  main  plasma.  The  implosion  of  the  main  plasma 
can  be  computed  with  the  resistive  magneto-hydrodynamic  equations,  with  proper 
includion  of  radiation  production  and  transport.  However,  in  the  computations 
the  exterior  of  the  plasma  sometimes4  behaves  in  an  unphvsical  way:  the  plasma 
temperature  outside  the  pinch  T  can  exceed  100  keV,  the  plasma  conductivity 
increases  excessively  when  using  the  standard  Spitzer  formula  a  oc  T3^2 ,  and  the 
high  conductivity  prevents  the  current  from  entering  the  z-pinch.  The  problem 
disappears  when  the  conductivity  is  calculated  properly,  keeping  into  account  the 
transition  from  a  plasma  dominated  by  collisions  to  a  collisionless  plasma  (Section 
4.3). 

In  this  strongly  magnetized,  collisionless  exterior  plasma  resistive  hydrody¬ 
namics  is  invalid.  Instead,  the  plasma  electrons  drift  in  the  confining  magnetic 
field  Bo(r )  and  the  combined  applied  and  space-charge  electric  field  E  =  (£b(r.r). 
Er(z,r)).  The  E:  x  Bo  drift  inward  generally  dominates.  The  current  density  j. 


no  longer  satisfies  Ohm's  law,  but  instead  is  found  by  setting  the  power  density 
j:E,  equal  to  the  change  in  kinetic  energy  density,  or  equivalently,  by  taking  the 
appropriate  moments  of  the  single  particle  distribution  function  The  various  cur¬ 
rent  losses  are  small  when  there  is  little  material  outside  the  z-pinch.  In  order  of 
magnitude,  Iioss  j  Itotai  ^  (Nou*  Tout  /A pinch.  "Epmch  )  ^  .  where  Vpmc/i.out  is  the  num¬ 
ber  of  particles  in  the  bulk  of  the  pinch  (outside),  and  Tplnch,jUt  's  the  applicable 
temperature. 

A  reasonable  mass  per  unit  length  for  the  z-pinch  bulk  is  perhaps  100 fig /cm. 
while  the  particle  density  in  the  z-pinch  periphery  n  ~  10 13 /cm3  (~  1  mTorr ):  over 
100cm2  the  mass  per  unit  length  in  the  outside  the  pinch  is  ~  0.01  fig /cm.  Setting 
the  typical  temperatures  about  equal  suggests  a  relative  current  loss  of  ~  1%,  see 
section  IV.  Propagation  of  the  electrical  pulse  through  the  exterior  plasma  isunder 
continuing  study  by  one  of  us  (RET). 

In  their  inward  motion  the  electrons  keep  the  angular  momentum  p$  constant, 
and  therefore  the  electron  velocity  v$  increases  as  1/r  and  the  parallel  pressure  p\\  as 
1/r2.  Also  the  magnetic  moment  p  =  mu2  /2 B  oc  pj_  / B  is  conserved.  Therefore  the 
perpendicular  pressure  pL  increases  with  the  magnetic  field  as  1/r.  The  difference 
between  the  parallel  and  perpendicular  pressure  increases  as  the  implosion  of  the 
collisionless  plasma  proceeds  toward  the  center. 

When  the  parallel  and  perpendicular  pressure  are  unequal  an  exact  self-similar 
oscillating  pinch5,8  is  no  longer  possible.  However,  the  oscillation  remaines  sepa¬ 
rable.  The  separable  pinch  has  a  low-density  domain  around  the  axis  in  response 
to  the  centrifugal  potential  that  represents  the  non- zero  angular  momentum  of  the 
electrons  (Section  4.5). 


4.2  ESTABLISHMENT  OF  MAGNETIC  INSULATION  AROUND  A  Z-PINCH. 

We  have  done  computations  on  a  generic  z-pinch  diode'  that  inspires  the  geo¬ 
metrical  and  pulse  power  parameters  in  the  present  work.  The  anode-cathode  gap. 

3  cm  in  this  case,  is  bridged  by  an  ideal  conductor  that  mimics  the  z-pinch.  The 
z-pinch  radius,  0.6  cm,  is  representative  of  multiple  wire  loads8  but  substantially 
smaller  than  the  initial  radius  of  a  gas  puff. 

In  the  computer  simulations,  performed  with  MAGIC9,  electrons  are  emitted 
from  all  cathode  surfaces  assuming  space-charge  limited  flow,  i.e.  emission  continues 
until  the  normal  electric  field  vanishes.  No  ion  emission  is  allowed. 

The  voltage  at  the  diode  entrance  rises  linearly  on  a  5  ns  timescale  to  over 
1  MV.  The  current  satisfies  1= V/ L  and  is  quadratic  in  time,  exceeding  100  kA  at 
the  end  of  the  simulation,  for  a  magnetic  field  at  the  z-pinch  edge  up  to  about 

4  Tesla  (40  kG). 
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The  current  in  free  electrons  starts  once  the  electric  field  in  the  diode  feed 
exceeds  200  kV/cm.  At  this  time,  t=0.5  ns,  the  current  through  the  z-pinch  is 
insufficient  for  magnetic  insulation,  and  the  free  electrons  simply  cross  the  diode  gap. 
Later  in  time  these  electrons  become  trapped,  and  evolve  until  t=4.4  ns.  the  end  of 
the  pulse,  the  electrons  are  distributed  in  space  as  in  in  vortices.  The  magnetically 
insulated  space  charge  sheath  in  the  current  feed  appears  to  stream  off  the  cathode 
into  the  diode  already  filled  with  previously  emitted  electrons.  The  sheath  is  clearly 
unstable,  rolling  up  into  vortex-like  structures  that  persist  as  they  move  toward 
the  center  conductor.  Qualitatively  similar  results  are  found  in  simulations  with 
different  parameters,  e.g.  a  slower  voltage  risetime,  or  a  smaller  z-pinch  inductance 
more  representative  of  a  gas  puff. 

There  is  an  exact  analogy  between  strongly  magnetized  non-neutral  electron 
beams  and  inviscid  fluid  flow10.  In  two-dimensional  Cartesian  geometry  with  con¬ 
stant  magnetic  field  Bz  the  charge  density  p/eoB,  corresponds  to  the  vorticitv 
G  =  V  x  v.  The  electron  fluid  drift  velocity  vd  can  be  derived  from  the  stream 
function  tv  as  Vp  =  ez  x  Vip,  where  V2tv  =  Qz- 

The  exact  correspondence  breaks  down  in  cylindrical  geometry  because  the 
magnetic  field  varies  with  radius  r,  Be  =  poI/2irr;  now 

„  /on  3(1  /Be) 

Vxv  =  ~p/e0Be  -  Er — ^ — .  (i 

The  last  term  is  small  compared  to  the  first,  except  perhaps  near  the  central  con¬ 
ductor  where  the  blob  size  is  comparable  to  the  geometrical  scale.  Hence  in  our 
diode  the  analogy  with  inviscid  shear  flow  is  largely  valid. 

The  magnetically  insulated  space  charge  sheath  coming  off  the  cathode  is  an 
unstable  vortex  sheath.  In  fluid  dynamics  it  is  well-known  how  the  instability 
develops,  namely  into  coherent  vortex  regions  reminiscent  of  the  space  charge  blobs 
seen  in  our  simulation.  The  blobs  are  stationary  entities  because  the  electrons 
remain  on  equipotential  lines.  The  drift  velocity  vD  =  ExB  is  perpendicular  to 
the  potential  gradient  E  =  -V<p,  and  vp  •  Vrt  =  0.  In  a  cylindrically  symmetric 
blob  the  charge  converted  by  the  drifting  electrons  is  then  invariant,  and  the  blob 
is  a  stationary  entity.  It  is  also  known  that  isolated  vortex  regions  are  relatively 
stable. 

Interaction  between  vortices  is  a  strong  mixing  process11.  The  counterpart  in 
our  simulations  should  be  a  uniform  spreading  of  the  space  charge  blobs  over  the 
vacuum  region  of  the  diode.  We  have  not  studied  vortex  mixing  because  the  pulse 
is  over  before  mixing  can  be  completed. 

The  leakage  current  in  free  electrons  peaks  at  about  2  ns  and  decreases  there¬ 
after.  The  leakage  current  should  be  well  approximated  by  the  drift  current  parallel 
to  the  cathode  in  the  thin  space  charge  sheath  that  rolls  off  the  cathode.  In  an 
insulating  magnetic  field  Be  the  sheath  electrons  drift  along  the  cathode  into  the 
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diode  with  the  drift  velocity  vz  ~  Er  x  Bg/Bg.  The  electric  field  varies  through¬ 
out  the  electron  sheath  from  zero  at  the  metal  surface  of  the  cathode  to  Er  at 
the  top  of  the  sheath,  but  the  magnetic  field  is  mainly  from  the  central  current, 
hence  Bg  ~  I/2~rt.  The  current  in  the  electron  sheath  7e  is  then  approximately 
/e  =  2; Trtcrv:.  The  surface  charge  density  a  of  the  cathode  sheath  is  cr  =  e0Er. 
where  Er  is  the  normal  electric  field  at  the  outside  of  the  electron  sheath.  The  sheath 
is  thin  compared  to  the  pulselength  width  d,  which  is  typically  small  compared  to 
the  inner  radius  r,  of  the  pulseline.  Then  the  electric  field  can  be  approximated  bv 
Er  =  V/d.  Usually  the  inductive  voltage  dominates  because  the  risetime  r  is  so 
fast,  and  V  =  LI  ~  LI/t  ( L  is  the  diode  inductance,  I  the  total  diode  current). 
With  these  approximations  the  relative  importance  of  the  electron  current  is 


where  Z0  =  (^0/eo)2  =  c  =  377  Q,  and  Zd  is  the  total  diode  impedance. 

The  literature  on  magnetic  insulation12-14  contains  estimates  along  the  same 
lines  as  the  above.  In  addition  it  is  shown  in  detail  that  the  single  particle  consid¬ 
erations  used  here  correspond  well  to  a  treatment  that  includes  the  self-electric  and 
-magnetic  fields. 

A  typical  z-pinch  diode  perhaps  3  cm  wide,  and  the  inductance  L  ss  6  nH .  A 
typical  pulse  risetime  may  be  r  «  12  ns,  for  a  diode  impedance  of  about  0.5  fi.  In 
contrast,  the  vacuum  impedance  Zo/2tt  =  60Q,  and  the  impedance  ratio  is  1/100. 
The  electron  current  /„  is  thus  always  a  small  fraction  of  the  total  current  /,  whence 
/e  can  be  ignored  for  imploding  z-pinches. 

For  single  exploding  wires,  in  contrast,  the  vacuum  current  may  become  com¬ 
parable  to  the  conduction  current.  As  an  example,  the  inductance  of  a  single  25/zm 
radius  wire  load  in  a  3  cm  wide  diode  is  about  45  nH.  For  a  fast  10  ns  risetime 
pulse  the  (actor  ( Z/?/Z0  ~  0.1.  but  a  relatively  small  gap  size  partly  compensates. 
When  r,jd  %  3,  lor  example,  /e / /  of  order  1/10.  In  this  case  there  should  be 
a  sizeable  leakage  current  in  free  electrons,  potentially  generating  a  measurable 
bremsstrahlung  x-ray  signal.  Such  a  signal  may  have  been  observed10  in  in  the 
Gamble  II  device  at  the  Naval  Research  Laboratory,  although  no  quantitative  mea¬ 
surements  are  available.  If  corroborated  by  further  work  the  bremsstrahlung  x-ray 
signal  may  be  an  excellent  indicator  for  the  arrival  time  of  the  pulse. 


4.3  TRANSITION  TO  COLLISIONLESS  PLASMA. 

The  electron  energy  distribution  function  in  the  periphery  of  a  radiating  z-pinch 
is  relatively  easy  to  calculate  by  assuming  that  the  electron-ion  collisions  isotropize 
the  plasma  electrons,  but  are  unimportant  in  energy  exchange  between  electrons. 


Then  the  isotropic  part  of  the  electron  distribution  function  /0  satisfies  a  highly 
simplified  form  of  the  Fokker-Planck  equation16,  viz.. 


Oja  _  eE  d  2 
dt  3m  t’2  dv  1 


C(foJo). 


(D1 


Here  C(fo-fo)  is  the  Fokker-Planck  collision  temr  for  multiple  scattering  between 
electrons,  / 1  is  the  anisotropic  part  of  the  electron  distribution  function  f  —  f0  - f 
fi  +  ■  ■  •  in  the  electric  field  direction  E/E, 


h  = 


rn^etU  + 


Here  is  the  electron-ion  collision  frequency,  and  the  magnetic  field  B  enters 
through  the  electron  cyclotron  frequency  ^ce  =  eB/mc.  The  perturbation  treat¬ 
ment  is  valid  provided  /i  remains  small  compared  to  /0. 


Quantitatively,  the  transition  to  a  collisionless  plasma  can  be  studied  with  a 
numerical  solution  of  the  isotropic  Fokker-Planck  equation,  implemented  according 
to  Ref.  IT.  The  distribution  function  changes  in  the  transition  from  a  collisional 
plasma,  with  low  temperature,  to  a  hotter  plasma  where  collisions  are  no  longer 
important.  Eventually  the  distribution  function  approaches  the  asymptotic  form18 
/(e)  ~  exp[— e5/2 /e^2  ( r )] :  for  large  e  the  distribution  function  /(e)  decreases  much 
faster  than  a  Maxwellian.  The  normalized  temperature  eo(r)  increases  with  time 
as  r2/5 ,  and  the  electrical  conductivity  decreases  with  time: 


The  temperature  and  conductivity  versus  time  eventually  reach  an  asymptotic  state. 

How  long  it  takes  the  plasma  in  the  z-pinch  periphery  to  reach  this  asymptotic 
state  depends  largely  on  the  initial  density.  To  our  knowledge  there  are  no  mea¬ 
surements  on  the  periphery  of  an  imploding  z-pinch.  The  initial  gas  density1'*  may 
be  on  the  order  of  n3^3  ~  1013  —  lO^/cm3.  100-1000  times  the  base  pressure.  If 
the  plasma  were  a  hypothetical  fully  ionized  neon  at  1  eV  the  collision  time  would 
be  1  ns.  and  the  electron  distribution  reaches  its  asymptotic  state  in  the  pulse  time 
(~  10ns).  For  fully  ionized  neon  at  100  eV  the  collision  time  is  1  /is.  whence  the 
distribution  function  can  not  change  during  the  pulse. 

If  a  peripheral  plasma  were  important  for  the  character  of  the  z-pinch  implosion 
there  should  be  some  unexplainable  difference  between  z-pinches  that  start  from  an 
injected  gas  on  the  one  hand,  and  multiple  wires  on  the  other.  Slight  differences  that 
have  been  observed  can  be  accounted  for  by  other  factors  than  a  current  leakage 
in  the  pinch  periphery.  For  example,  each  of  the  wires  in  a  multiple-wire  implosion 
blows  up  early  in  the  pulse.  Sometimes20  the  blow-off  plasma  is  torn  off  and  im¬ 
plodes  on  axis  even  before  the  bulk  of  the  original  wires  has  moved.  The  lack  of 
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positive  experimental  evidence  suggests  but  does  not  exclude  that  a  plasma  outside 
an  apparently  well-defined  current  channel  has  a  negligible  effect  on  the  pinch. 

The  remainder  of  this  paper  assumes  that  the  z-pinch  periphery  contains  a  col¬ 
lisionless  plasma,  with  density  and  temperature  to  be  determined.  The  magnitude 
of  these  parameters  will  determine  whether  the  peripheral  plasma  is  important  in 
z-pinch  implosions. 


4.4  COLLISIONLESS  CURRENT  DENSITY. 

In  a  collisional  plasma  the  collisions  are  the  dominant  influence  that  balances 
the  motion  of  the  plasma  electrons  against  the  acceleration  by  the  elect  ric  field.  The 
resulting  collisional  electron  drift,  and  thus  the  current  density  j,  is  proportional 
to  the  electric  field  E.  namely  j  =  <rE.  The  electrical  power  input  into  the  plasma 
ends  up  in  increased  random  velocities  of  the  plasma  particles,  i.e..  Ohmic  heating. 

Even  in  the  absence  of  collisions  the  magnetized  plasma  in  the  z-pinch  periphery 
can  absorb  energy  from  the  electrical  pulse.  In  this  case  the  energy  does  not  go  into 
random  motion  of  the  plasma  particles,  but  instead  is  put  into  organized  motion. 
This  motion  was  identified  but  not  completely  amalyzed  under  the  present  contract. 
However,  much  subsequent  work  has  shown  some  interesting  phenomena.  Some  of 
this  is  included  here  for  completeness,  and  to  give  a  perspective  on  how  research 
changes  direction  in  response  to  unanticipated  results. 

A  strongly  magnetized  plasma  in  an  electric  field  is  not  accelerated  by  the  field, 
but  drifts  with  a  velocity  v  =  E/ B.  For  constant  electric  field  the  energy  remains 
constant,  but  when  the  electric  field  increases  in  time  the  plasma  accelerates.  The 
kinetic  energy  density  pv2  increases  also  ( p  is  the  plasma’s  mass  density),  and 
electrical  power  is  absorbed.  Likewise,  electrons  that  E  x  B  drift  into  a  region  of 
increased  magnetic  field  spin  up.  increasing  their  perpendicular  energy  u.’j_  in  an 
effort  to  keep  t  he  magnetic  moment  p  ~  w<_/B  constant.  The  increase  in  u'<_ .  or 
the  increased  magnetization,  is  another  power  sink. 

The  angular  momentum  about  the  axis.  P,)  =  mri'g.  is  conserved  in  a  cylindri- 
cally  symmetric  geometry,  and  U||  increases  as  1/r.  A  drift  toward  the  axis  increases 
the  kinetic  energy  parallel  to  the  magnetic  field  Wjj  =  mvjj/ 2  as  1/r2. 

Energy  conservation  is  then  expresed  by  the  formula 

j-E  =  d(I<E)/dt  +  v-F.  (12 

where  the  pressure  force  is  F  =  — V.p.  and  p  is  the  pressure  tensor. 

The  macroscopic  quantities  that  enter  into  the  fluid  equations  are  averages 
over  the  particle  distribution  function.  Without  collisions  the  distributions  are  not 
necessarily  maxwellian,  but  the  exact  functional  form  is  unimportant  because  only 
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the  averages  count.  Instead  of  the  distribution  in  terms  of  the  velocity  components 
in  cylindrical  geometry  it  is  better  to  use  the  distribution  function  in  terms  of  the 
dynamical  invariants,  the  magnetic  moment  p  and  the  angular  momentum  Pg.  The 
distribution  function  is 


/(r,  v,  t.) 


m  2  r  "] 

■  2x  B(r.t)  J 


n(r,  t)  F(fi,  Pg.  r) 


i  i  ■?  \ 

\  A  ’J  I 


The  factor  in  [•  •  •]  is  the  Jacobian  of  the  transformation  from  the  velocity  vector 
v  to  the  gyrophase  angle  <5,  magnetic  moment  p.  and  angular  momentum  about 
the  axis  Pg.  In  the  absence  of  collisions  the  magnetic  moment  is  approximately 
conserved,  but  even  with  collisions  Pg  is  invariant.  Then  the  distribution  function 
F  is  constant  in  time.  The  particle  density  n(r,t)  is  no  longer  part  of  F.  which  is 
normalised  as  _ 

mQ,  r  TT  1  /  dpdPg  F(p.Pg.r)  =  1.  (14) 

27T  B(r,  t)  J  J 

for  both  electrons  and  ions.  With  the  assumption  of  quasineutrality  the  mass  density 
P  is  P  =  {rrielectrona  +  ml0n,  )n(r,t). 

The  parallel  pressure  is 

PH  =  pL\  =  J  dpdPg  FPg  /m2,  (15) 

and  the  perpendicular  pressure  becomes 

Pl  =  pWi.  =  J  dpdPg  Fp/m ,  (16) 


and  the  summation  over  species  is  understood. 

The  axial  current  density,  calculated  either  by  momentum  balance  or  by  aver¬ 
aging  of  ntv,.  is  then 


pc2  D(E/B)  c  pL\\  I  (){ rpB\VL  ) 
B  Dt  B  .  r3  r  dr 


(IT) 


I  he  first  term  is  related  to  the  acceleration  of  the  E  x  B  drift  motion,  the  second 
comes  from  angular  momentum  conservation,  or  /?||,  and  the  third  is  the  magne¬ 
tization  contribution,  from  pL.  In  the  acceleration  term  the  ions  dominate:  the 
electrons  are  more  important  for  the  two  pressure-derived  terms  without  thermal 
equilibrium  between  electrons  and  ions.  The  effect  of  Eq.  (17)  on  electromagnetic 
energy  transport  is  presently  under  study  by  one  of  us  (RET). 
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4.5  EXACT  SEPARABLE  IMPLOSION  WITH  pL  ±  P]]. 


For  special  initial  conditions  there  exists5,6  an  exact  solution  to  the  fluid 
equations  that  describe  strongly  magnetized  plasma  in  the  periphery  of  a  current- 
carrying  z-pinch.  These  equations  are  standard. 


dn/dt  +  V  ■  (nv)  =  0. 


i  18) 


and 


v  =  cE  x  B /B2, 
dB/dt  =  V  x  (v  x  B), 


V  x  B  =  4jrj/c. 


(19) 

(20) 


(21 


However,  in  what  follows  the  momentum  equation  includes  an  anisotropic  pressure 
*pt.  pdv/dt  —  — V  •  +  j  x  B/c. 

An  isotropic  pressure  implies  that  the  collisions  are  sufficient  to  keep  the  veloc¬ 
ity  components  randomized.  This,  then,  implies  a  finite  conductivity,  in  contrast  to 
the  assumption  of  infinite  conductivity  which  demands  a  collisionless  plasma.  How¬ 
ever,  as  the  plasma  gets  closer  to  the  axis  in  a  cylindrical  geometry  the  pressure  pq 
parallel  t  the  magnetic  field  B  =  Bg{r)  increases  as  1/r2,  while  the  perpendicular 
pressure  p\_  increases  as  1/r.  In  fact,  in  any  curvilinear  magnetic  field  an  isotropic 
pressure  is  strictly  speaking  inconsistent  with  infinite  conductivity. 

The  fluid  equation  of  motion  in  the  radial  direction,  for  an  anisotropic  pres¬ 
sure,  includes  an  centrifugal  force  /r3  that  reflects  the  conservation  of  angular 
momentum, 
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The  centrifugal  term  is  L\ i  =  L\\eiectTona 
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r  dr 
+  £|i,,ons  •  and  for  each  species 


(22) 


pL^ j  =  m  J  d3vfr2u 
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where  ug  is  the  particle  velocity  along  the  magnetic  field.  In  terms  of  manifestly 
invariant  quantities 

pL\  =p  J  dpdPg  F-^2-  (236) 

Likewise,  the  perpendicular  pressure  is  connected  to  the  kinetic  energy  perpen¬ 
dicular  to  the  magnetic  field. 


pH  L  =  pBW L 


(24n) 
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which  is,  in  terms  of  the  invariants. 


pBWL  —  pB  J  (lj.td.Pg  Fjl/r 


Equation  (22)  allows  a  selfsimilar  oscillation  of  the  pinch  plasma3  when  the 
angular  momentum  L\\  vanishes,  but  even  for  non-zero  L\\  a  self-similar  oscillation 
is  possible.  A  comprehensive  discussion  of  these  solutions  is  deferred  to  a  future 
paper21 ,  but  here  a  simple  example  is  shown. 

Henceforth  the  dependent  variables  B.  p ,  L p ,  W L  and  their  derivatives  are 
functions  of  the  single  selfsimilar  variable  £  =  r/a(f),  where  a(t)  is  a  typical  scale 
length  that  contains  the  sole  time  dependence.  Then  Eq.  (22)  becomes 

d2a  '*,’m  u,'2  uj2, 

TT"  =  4  +  (25] 


where  a(t)  is  the  normalized  scale  length  a(t)  —  a(t)/a0.  Non-vanishing  angular 
momentum  is  responsible  for  the  term  with  1/a3,  with  u;||2(£)  =  £||(0/£4- 

In  general  the  coefficients  uj2  =  u2(£)  depend  on  the  spatial  variable  f,  and  Eq. 
(25)  does  not  separate.  Separation  is  feasible  only  when  u>jj  is  constant.  Recall  that 
L2 j  is  an  average  of  where  Pg  ~  vsr  is  the  particle  angular  momentum.  Since 
vg  ~  u tgr  also  Pg  ~  u jgr2 ■.  constant  u>^  then  implies  a  parallel  temperature  propor¬ 
tional  to  (radius)2,  and  an  initial  angular  momentum  proportional  to  (radius)4. 

The  perpendicular  pressure  term  gives  rise  to 


U7,  = 


1  d[HL(0(p] 


which  should  no  longer  be  a  function  of  £.  In  this  term  the  spatial  dependence 
of  the  density  p(£)  can  be  played  olf  against  the  perpendicular  thermal  energy 
H l  (t)  =  (£)•  Either  p  or  HL  can  be  chosen  freely.  It  is  easiest  to  separate 

the  spatial  dependence  in  H L  by  setting  Hx(t)  =  H0/g{£).  The  mass  density  is 
then 


t  l 


The  acceleration  term  contains  the  magnetic  field  through  u:\  =  B2 /\irpa^. 
vhich  is  the  transit  frequency  of  an  Alfven  wave  though  the  scale  size  a0.  A  constant 
therefore  implies  a  specific  magnetic  field  given  by 


d(ZB)2 
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where  B p0  is  the  magnetic  field  at  the  characteristic  distance  a0.  Equa¬ 
tion  (28)  is  really  an  expression  for  the  total  current  enclosed  by  the  pinch  within 
a  distance  because  £B  =  2/(£)/ca0. 

The  function  g(£)  is  arbitrary  provided21  that  the  density  p  be  positive  and 
finite  everywhere,  and  that  the  various  u.,2's  remain  positive.  The  simplest  function 
to  satisfy  these  requirements  is  g(£)  =  £2.  Then  the  various  integrations  can  be 
done  in  closed  form.  For  example,  the  density  becomes 

p  =  potexp  -  (£4  -  1  )/4,  (29a) 

where  p0  is  given  by  the  mass  per  unit  length  within  £.  Then  the  mass  per  unit 
length  within  £  is  given  by 

2 -alj^dUP  =  N/2Ta2poe1/47(-l/4.^/-i)-  (296) 


where  7 (a,x)  is  the  incomplete  gamma-function.  Likewise,  the  current  carried  by 
the  corona  out  to  £  becomes 


I2(0-I\ 0) 


c2°-lBlo 


>/2 


7(1/4, £4/4), 


(30) 


where  7(0)  represents  the  current  carried  by  a  collisional  pinch  located  on  axis.  This 
is  reminiscent  of  the  Bennett  pinch  relation,  72  a  NT,  for  good  reason:  The  corona 
contributes  its  pressure  to  the  current  carrying  capacity  of  the  pinch. 


4.6  SUMMARY  AND  CONCLUSIONS. 

Current  can  be  carried  in  the  periphery  of  a  z-pinch  by  various  mechanisms. 
Initially  the  z-pinch  carries  no  current,  and  the  inductive  voltage  at  the  edge  of  the 
cathode  awav  from  the  pinch  material  can  drive  a  currrent  in  free  electrons.  The 
free  electrons  can  be  the  bulk  of  the  current  at  the  very  start,  but  this  current  is 
cut  off  rapidly  once  the  current  starts  flowing  through  the  z-pinch. 

The  tenuous  plasma  outside  the  z-pinch  can  also  carry  current.  As  a  remnant  of 
background  gas  or  perhaps  the  edge  of  the  injected  gas  column  in  a  puff-gas  z-pinch. 
this  plasma  is  not  well-specified  either  experimentally  or  theoretically.  How  much 
current  is  shunted  through  this  plasma  depends  linearly  on  the  plasma  density,  and 
also  on  parameters  such  as  the  temperature  and  the  magnetic  field  in  the  periphery, 
i.e.  the  current  in  the  main  z-piuch. 

The  plama  temperature  in  the  periphery  is  computed  from  a  simplified  Fokker- 
Planck  analysis,  assuming  that  the  electron  velocity  remains  isotropic  through  col¬ 
lisions  with  the  ions.  Plasma  heating  proceeds  until  the  electron  collision  time  is 
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comparable  to  the  pulse  time,  resulting  in  a  temperature  limit  of  something  like  50 
eV  (weakly  dependent  on  plasma  density). 

The  pressure  of  the  collisionless  plasma  in  the  periphery  becomes  aniso.  rooic 
as  the  plasma  contracts  toward  the  axis,  due  to  the  separate  conservation  of  the 
electron's  magnetic  moment,  and  the  electron's  angular  momentum  about  the  axis. 
Under  these  conditions  the  plasma  can  still  oscillate  in  a  self-simiiar  manner,  and 
for  constant  z-pinch  current  and  balanced  oscillation  frequencies  (uqj .  ^  )  the 

periphery  can  be  in  an  exact  equilibrium. 

As  the  current  increases  in  the  power  pulse  part  of  its  energy  is  siphoned  off  by 
the  periphery.  How  much  depends  on  the  plasma  parameters,  which  are  unknown. 
For  a  tenuous  periphery,  n  <  10 13 /cm3,  with  the  corresponding  temperature,  T  ~ 
o.OeV’,  little  power  is  diverted  from  the  z-pinch.  For  densities  appropriate  to  edge 
gas,  n  ~  1 0 1 5 / cm3 ,  or  higher  temperatures  T  ~  oOeU.the  power  loss  is  substantial. 
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SECTION  5 


LOW-INDUCTANCE  PLASMA  RADIATION  SOURCE  LOADS 

Compact  high-power  electrical  systems  based  on  inductive  energy  storage  for 
driving  plasma  radiation  source  loads  depend  on  fast  opening  switches  that  must 
break  many  megamperes.  This  is  no  simple  task  even  in  the  best  of  circumstances, 
but  it  is  especially  difficult  when  opening  the  switch  generates  a  high  voltage.  An 
ideal  opening  switch  would  be  able  to  decrease  the  current  I  to  zero  in  a  time  r. 
irrespective  of  the  voltage  this  process  generates.  When  the  current  is  commuted 
into  an  inductive  load  L  there  appears  a  voltage  spike  with  maximum  Vm  LI jr 
(in  addition  to  any  other  voltage,  e.g.  a  resistive  voltage  V  =  RI  ).  In  other  words, 
when  the  current  is  not  immediately  diverted  into  a  load  the  current  will  continue 
to  flow  in  the  switch.  Depending  on  the  physical  characteristics  of  the  switch  the 
consequent  energy  deposition  will  modify  the  switch  operation,  usually  for  the  worse. 
For  example,  the  switch  may  restrike,  or  not  open  at  all. 

The  punishment  the  opening  switch  is  getting  can  be  reduced  by  minimizing 
the  inductance  of  the  load  circuit.  Part  of  the  minimisation  involves  simply  geom¬ 
etry.  A  physical  separation  between  the  switch  and  the  load  implies  an  inductance 
L  %  g0i  x  G,  where  i  is  a  typical  length  scale  such  as  the  switch-load  distance, 
and  G  is  a  geometrical  factor.  Clearly,  it  is  advantageous  to  make  the  switch-load 
distance  as  small  as  possible.  One  way  to  do  this  is  by  putting  additional  switches 
in  the  system  that  divide  the  inductance  among  themselves.  Here  we  describe  some 
ideas  that  may  be  useful  in  this  general  context. 

A  low-inductance  wire  load  can  be  made  by  stretching  wires  between  the  outer 
edge  of  anode  and  cathode,  followed  by  rotating  one  of  the  electrodes  over  an  angle 
O.  No  twist,  or  0  =  0,  gives  a  standard  wire  load,  and  o  <  v  creates  the 
so-called  X- pinch  recently  investigated  by  Soviet  workers.  Any  angle  in  between 
0  and  x  leads  to  wires  with  a  hyperbolic  envelope.  Hence  there  is  a  continuous 
transition  between  the  standard  configuration  and  the  X-pinch.  The  envelope  can 
be  expressed  in  terms  of  the  twist  angle  o.  but  the  simplest  way  of  doing  this  is  in 
terms  of  the  waist  radius  r0  and  the  asymptotes  with  angle  tg  0  —  r0/'2/?0,  where 
r0  is  the  diode  gap  length  and  RQ  the  electrode  radius 

It  is  clear  that  the  twisted  wires  have  less  inductance  than  wires  strung  straight 
between  the  electrodes  at  the  radius  r0  of  the  hyperboloid’s  waist,  but  a  larger 
inductance  than  wires  at  the  electrode  edges  Rq .  Optimum  radiation  production 
may  require  that  the  mass  per  unit  length  fi  —  m/i  and  the  radius  r  are  related  by 
gr2  —  C.  where  C  is  a  (machine-dependent)  constant.  In  this  case  the  inductance 
of  the  twisted  wires  is  undoubtedly  more  favorable  than  the  straight  wires. 

Twisting  the  wires  only  in  the  positive  direction  will  generate  a  magnetic  field 
Bz  along  the  machine  axis:  this  can  be  cancelled  by  twisting  the  same  number  of 
wires  in  the  negative  direction. 
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As  usual  in  pulse  power  it  seems  much  simpler  to  try  out  the  geometry  than 
to  simulate  it  numerically,  especially  if  the  angle-dependent  features  are  to  he  kept. 
Assuming  that  the  wires  behave  as  their  hyperboloidal  envelopes  makes  life  suffi¬ 
ciently  simple  to  allow  simulation  on  a  2-D  hydromagnetics  code. 

Since  the  load  geometry  is  a  continuous  function  of  0  it  may  be  possible  that 
the  radiation  is  maximized  at  a  value  intermediate  between  0  and  ~.  e.g.  at  o0- 
Then  the  load  concept  has  merit. 

As  with  straight  wire  loads  the  radiation  yield  may  be  further  increased  by 
letting  the  mass  per  unit  length  change  with  axial  position,  i.e  /i  =  with 

larger  mass  in  the  pinch  center.  However,  the  twisted  wires  load  may  be  particularly 
suited  for  use  in  generators  based  on  inductive  energy  storage,  where  the  load 
inductance  may  make  the  difference  between  a  failing  or  a  succesful  switch. 

Another  way  to  decrease  the  switch-load  distance  is  by  putting  an  intermediate 
switch  in  the  space  between  the  two.  This  can  be  done  with  wire  loads.  The 
standard  PRS  load  is  halfway  down  toward  the  center,  but  there  is  a  second  wire 
load  on  the  outside  of  the  electrodes  where  the  switching  takes  place,  at  Ro.  These 
wires  should  have  the  appropriate  mass  per  unit  length  /i  for  them  to  reach  the 
radial  position  of  the  PRS  load  wires  at  the  desired  switching  time.  Presumably, 
the  outer  ’’switch”  wires  divert  the  current  from  the  load  wires  because  of  their 
smaller  inductance.  In  the  limit  of  infinitely  many  wires 

L  switch.  1  rt(Rsujitch  /Ro)  ^ 

Lprs  In  (Rprs /Ro) 

which  is  approximately  A  R/ Ro  1:  for  a  finite  number  of  wires  the  current  is 

not  as  well  diverted.  In  this  more  realistic  case  there  is  leakage  of  the  magnetic  field 
through  the  gaps  between  the  switch  wires,  which  adds  to  the  inductance. 

The  material  from  the  switch  wires  will  not  disappear  but  instead  it  may  collide 
with  the  PRS  wire  load,  or  it  may  gather  and  cook  on  axis.  To  avoid  any  interference 
in  radiation  production  by  the  switch  material  its  Z  should  be  small  compared  to 
the  Z  of  the  PRS  load.  e.g.  2,wltch  <  6:  plastic  or  graphite  should  work. 

For  a  change  it  seems  easier  to  evaluate  a  moderately  realistic  simple  0-D  model 
of  such  a  switch  enhancement  than  to  do  the  experiment,  but  we  haven't  done  it  as 
yet.  First  we  want  to  see  what  other  people  think  of  it. 

The  (solid  or  liquid)  foil  switch  may  become  an  interesting  and  useful  topic  for 
numerical  modeling.  In  this  switch  a  foil  is  blown  across  a  small  cathode-anode  sap 
by  j  x  B  forces.  In  the  models  we  are  familiar  with  the  bulk  of  the  switch  foil  is 
presumed  to  remain  solid,  or  at  most  it  may  become  liquid,  but  there  are  little  or  no 
plasma  effects.  However,  at  the  edge  of  the  foil  where  it  touches  the  inner  (cathode) 
part  of  the  pulseline  it  is  hard  to  see  how  the  foil  can  avoid  becoming  a  plasma: 
especially  when  the  foil  starts  to  shoot  over  the  edge  in  its  opening  act  there  are 


mega-amperes  flowing  in  an  ever-thinning  sheath  of  material,  which  is  bound  to 
become  a  highly  conducting  plasma  that  streaks  across  and  into  the  cathode-anode 
gap.  exactly  what  we  are  trying  to  avoid. 

It  seems  crucial  to  switch  performance  how  much  plasma  gets  into  the  gap:  if 
there  is  little  the  plasma  should  collapse  on  axis  similar  to  our  switch  wires,  but  if 
there  is  too  much  an  arc  can  be  struck  between  the  electrodes,  and  the  switch  never 
opens,  i.e.  never  transfers  the  current  to  the  load.  We  propose  to  study  this  part 
of  the  process  more  in  depth,  partly  by  using  existing  hydrodynamic  codes  with 
addition  of  the  appropriate  conductivity  expressions.  Another  part  of  our  efforts 
will  be  to  use  literature  data  for  help  in  identifying  the  processes  that  dominate  the 
arc  plasma. 

Past  experience  with  arcs  suggests  that  a  comprehensive  modeling  of  all  the 
events  in  the  "opening”  cf  the  foil  switch  will  be  exceedingly  difficult,  but  it  is  still 
necessary  to  attempt  exactly  this. 
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SECTION  6 


ELECTRON  ENERGY  DISTRIBUTION  IN  A  Z-PINCH 

6.1  INTRODUCTION. 

State  of  the  art  computations  on  the  radiation  production  of  imploding  z- 
pinches  follow  the  standard  hydromagnetic  laws  for  conservation  of  mass  and  mo¬ 
mentum,  augmented  by  a  complete  accounting  of  the  energy  density  as  this  changes 
due  to  compression.  Joule  heating,  and  radiation  cooling.  The  plasma  parameters 
are  typically  such  that  particle  collision  times  are  much  shorter  than  macroscopic 
timescales,  and  on  this  basis  the  electron  distribution  function  f(v)  is  assumed  to 
be  in  equilibrium,  i.e.  f(v)dv=  4xf(v)v2dv  is  isotropic  in  direction  and  Maxwell- 
Boltzmann  in  energy.  Electric  fields  or  spatial  gradients  usually  produce  small 
deviations  from  Maxwellian,  and  the  corresponding  coefficients  can  be  taken  from 
the  literature22,23 .  Also  the  atomic  excitation  rate  are  computed  by  integrating  the 
excitation  cross  section  with  a  Maxwellian  electron  distribution. 

In  the  plasma  periphery  where  the  plasma  density  is  thought  to  be  low  (ne 
~  1012-1015/cm3)  and  the  electron  temperature  to  be  high  (Te  1  keV  ?)  the 
collision  time  is  large  compared  to  the  pinch  time  (  100  /xs  vs.  100  ns).  Clearly 
the  electron  distribution  is  not  necessarily  in  equilibrium:  in  fact,  f(v)  behaves  as 
exp(—  v5/5v*h),  and  the  plasma  temperature  v2^  increases  in  time  as  t2/,s.  The 
electrical  conductivity,  defined  by  the  temperature  rise  divided  by  the  electric  field 
energy  E2,  then  decreases  with  temperature  as  T~3^2 .  This  result  solves  a  difficulty 
that  occurs  occasionally  in  the  pinch  periphery  using  the  radiation-MHD  equations. 

In  the  edge  of  the  plasma  the  density  appears  sufficient  to  suspect  a  Maxwell- 
Boltzmann  equilibrium  distribution,  but  inelastic  collisions  can  give  significant  de¬ 
viations  from  a  Maxwellian.  Inelastic  collisions  are  typically  less  frequent  than 
electron-electron  collisions,  v.  <C  uee.  However,  transitions  between  the  first  ex¬ 
cited  state  and  the  heliumlike  ground  state  in  moderate  atomic  number  elements 
(from  neon  to  iron  and  beyond)  produce  the  most  interesting  radiation.  Helium-like 
ions  can  typically  not  be  excited  by  the  thermal  electrons  in  the  bulk  of  the  dis¬ 
tribution,  but  only  by  the  electrons  in  the  distribution’s  tail.  Elastic  collisions  are 
relatively  infrequent  there  and  thus  ionization  has  a  surprisingly  large  and  depress¬ 
ing  effect  on  the  distribution  function’s  tail.  Consequently,  the  actual  excitation 
rate  may  be  lower  than  calculated  from  a  Maxwellian2"4,  but  the  amount  of  change 
depends  on  how  the  distribution  functions  are  compared. 

When  the  electron  distribution  function  is  allowed  to  change  its  temperature 
in  response  to  the  inelastic  collisions  the  inelastic  reaction  rate  changes  with  an 
appreciable  factor  that  reflects  the  reduction  in  plasma  temperature.  However,  with 
an  external  electric  field  to  keep  the  temperature  constant  the  reaction  rate  is  almost 
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the  same  as  without  inelastic  collisions.  In  both  cases  the  distribution  function's 
tail  beyond  the  inelastic  threshold  is  significantly  depressed  (the  computations  up 
to  now  use  a  single  inelastic  process). 

Smaller  numbers  of  electrons  at  energies  higher  than  the  main  inelastic  thresh¬ 
old  should  be  observable  by  reduced  radiation  from  more  energetic  excited  states 
compared  to  the  lower-lying  states,  and  there  may  be  inconsistencies  between  the 
plasma  temperature  determined  from  the  different  states  on  the  basis  of  radiative 
rates  computed  from  a  maxwellian. 

It  should  be  emphasized  here  that  the  pinch  axis  is  explicitly  excluded  from  the 
considerations  for  the  moment.  Along  the  axis  the  magnetic  field  does  not  prevent 
electrons  from  running  away  in  any  strong  electric  field  that  may  exist.  Electron 
runaway  will  be  studied  with  a  Monte  Carlo  computer  code. 


6.2  FORMALISM. 


The  plasma  of  interest  is  spatially  homogeneous  and  highly  ionized,  with  an  av¬ 
erage  ionization  state  Z=10  or  higher.  Electron-ion  collisions  are  perhaps  10  x  more 
frequent  than  electron-electron  collisions.  Energy  exchange  between  electrons  pro¬ 
duces  the  equilibrium  Maxwellian  shape  of  the  electron  distribution  function  after 
roughly  5  self-collision  times25 .  In  contrast,  isotropy  in  the  electron  velocity  vector 
is  governed  by  elastic  electron-ion  collisions,  and  the  distribution  function  should 
be  almost  isotropic  while  it  relaxes  toward  equilibrium.  Electron  fluxes  come  from 
small  deviations  in  isotropy:  because  of  the  fast  electron-ion  rate  any  anisotropy 
is  assumed  to  be  in  equilibrium  with  its  driver.  Hence  the  fastest  timescale  is  the 
electron- ion  collision  time  ret-  =  ,  while  the  timescale  for  changes  in  plasma 

temperature  may  be  on  the  order  of  the  electron-electron  collision  time  ree.  Unfor¬ 
tunately,  the  collision  times  are  not  very  different.  re,  ~  ree/Z,  and  the  region  of 
validity  of  the  approach  taken  here  may  be  more  limited  than  desirable. 

The  purpose  is  to  define  a  reasonably  complete  but  still  analytically  tractable 
model  for  the  the  electron  energy  distribution  function  in  the  presence  of  an  electric 
field  E.  magnetic  field  B.  and  inelastic  collisions.  Simplifications  include  the  neglect 
of  any  spatial  gradients,  the  effect  of  magnetic  fields  on  the  collision  integrals,  energy 
exchange  in  elastic  collisions  between  electrons  and  ions,  the  absence  of  plasma 
instabilities,  etc.  The  plasma  ionization  rate  and  electrical  conductivity  can  then 
be  estimated,  and  modifications  herein  may  motivate  further  numerical  work.  Table 
l  contains  the  some  relevant  parameters  for  the  plasmas  of  interest. 

The  equation  for  the  anisotropic  part  of  the  distribution  function  /i(v)  is2'1'26'2' 


eE  df0 
me  t/c,  dv 


(32) 
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where  i/ei  is  the  electron-ion  collision  frequency, 

nZY 
=  — 5~ 


(33) 


The  known  symbols  have  their  usual  meaning,  and  Y  measures  the  effective  ionic 
size  in  multiple  scattering  with  the  Coulomb  interaction,  Y=Air(ei /Ant^nif.  )2(n\. 
with  A  ~  9np  and  no  the  number  of  electrons  in  the  Debye  sphere.  Typically. 
fnA  ~  10.  In  terms  of  fundamental  quantities  Y=  c 4  (&irr2/?>)  fnA/6,  with  re 
is  the  classical  radius  of  an  electron  and  the  term  in  parentheses  is  the  Thomson 
cross  section.  Electron-electron  collisions  can  be  included  properly  at  tremendous 
expense  in  complexity  and  approximately  by  replacing  Z  by  Z+l.  Neither  correction 
is  necessary  since  electron-ion  collisions  dominate.  The  even  less  frequent  inelastic 
collisions  are  also  ignored. 

The  magnetic  field  B,  always  taken  perpendicular  to  the  electric  field  E.  can 
be  included  easily  by  introducing  an  effective  electric  field  E ejj  =  E/(l  -f  / vf ). 
with  u)ce  =  eB/m.  The  anisotropic  part  of  the  electron  distribution  function  is,  in 
these  approximations,  determined  explicitly  from  the  isotropic  distribution  function 
fo- 

Relaxation  toward  equilibrium  of  fo  follows  the  equation  (e.g.  Ref.  27,  eq. 
7.71b,  or  Ref.  24,  5.16) 


dfo 

dt 


eE  6 
3 mev2  dv 


{v2  f\ )  =  So  + 


Si, 


(34) 


where  the  electron-electron  energy  exchange  is  contained  in  So,  and  the  inelastic 
electron-ion  collisions  are  represented  by  S,.  So  is  given  by 


**Yn  d ,  r  1  r  2  2  u2  p  /„  2  df0] 

~  T"l(  /  fo  1  ")  +  f  ~r~  /  (!  fo  t’  dv  +  —  /  —  r  dv)— —  . 

Ov  J  0  -3c  T)  3  J v  v  dr 

( 35u ) 

Only  inelastic  collisions  with  large  energy  loss  are  considered,  and  therefore  S, 
is  simply 

5,  =  —i/,  fo  -f  replacement  term ,  (35 b) 


where  the  electron  replacement  term  is  such  that  the  number  of  electrons  is  con¬ 
served.  The  ionization  frequency  vt=nmcr,v,  with  n’  the  density  of  the  ions  that  can 
be  excited,  and  <r,  the  velocity-dependent  excitation  cross  section. 


Normalization  results  in  the  equation 


r)h 

dr 


fffCh  +  (  D  + 


e2  c  dh , 

e2  1  +  dC 


•s (<f)/i  +  replacement  term. 

(36) 
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Here25  £=v/v0  and  Vg=3kT/m.  The  distribution  function  h(£,  r)=fo(v,t)v'o. 
and  time  is  scaled  as  r=  t/t0,  where  t0=VQZ/nY=Z/ree.  basically  the  electron- 
electron  collision  time  at  the  thermal  speed  v0.  The  electric  field  E  is  normalized 
by  Eo,  which  is  the  field  that  in  a  collision  time  t0  produces  a  velocity  change  on  the 
order  of  vq,  eEoto/m  —  v q  3 Z .  The  normalized  cyclotron  frequency  0  —  ^ce  to/Z 
and  typically  Q  1. 

Relaxation  due  to  electron-electron  collisions25  is  described  by  the  C(£)  and 
D(£)  terms  only,  with 

=  f  drjTj2h .  (37) 

Jo 

1  /■«  £2  [°° 

D  =  —  /  dr/  r)4h  +  I  di]  T)h.  (38) 

Numerical  solution  indicates25  that  an  arbitrary  initial  distribution  function  be¬ 
comes  Maxwellian  in  approximately  5  Spitzer  collision  times. 

The  replacement  term  will  not  be  used  in  this  form  of  the  equation,  but  will 
enter  to  conserve  electron  number  when  seeking  a  stationary  distribution  function. 

An  alternate  form  of  the  same  equation  is  obtained  by  using  the  (normalized) 
energy  e  =  £2, 

^57  =ah(£)/  +  /?(£)§71  -  Q(e)/  +  “a/a-  m 

Here  the  coefficients  are 


II 

Q 

II 

MJ 

/V/*, 

(40) 

■)  =  f  ^  /<*  +  ^  I 

The  ionization  loss  coefficient  a  is 

r  oo 

f  dt  + 

E 2  e3 

(41) 

e-  i  +  n  -E ' 

n-  T  4Qu(e) 
n  Z2tR  £nA 


where  is  the  collision  strength  for  the  atomic  levels  of  interest,  and  cr  is  the 
Rydberg  energy,  13.6  eV.  The  collison  strength  is  connected  to  the  inelastic  cross 
section  bv 


<r" 


7rar 


a, 

z2 


( R 
t 


(43) 


and  any  multiplicity  of  levels  is  suppressed.  The  excited  state  ionic  density  is  n*.  In 
this  form  of  the  equation  the  replacement  coefficient  becomes  simply  Qjk(e)  =  a(e  + 
t,-),  where  c,  is  the  excitation  energy. 


Eventually  the  single  inelastic  cross  section  will  be  replaced  by  a  summation 
over  all  the  processes  going  on,  including  the  super-elastic  collisions  that  put  en¬ 
ergetic  electrons  back  into  the  plasma,  and  collisions  that  absorb  little  energy  per 
interaction  but  that  collide  fairly  often.  For  the  present  purpose  these  processes  are 
distracting,  and  they  will  be  dealt  with  later. 

The  following  treats  two  special  cases  of  Eq.  (39),  the  first  relating  to  the 
plasma  periphery  where  the  collision  energy  exchange  terms  C  and  D  are  appropriate 
(though  negligible),  and  the  distribution  function  changes  in  time.  The  second  case 
treats  the  plasma  in  the  pinch  edge  away  from  the  axis,  where  the  distribution 
function  reflects  a  stationary  balance  between  collisional  energy  redistribution  and 
ionization. 


6.3  PINCH  PERIPHERY. 

In  the  pinch  periphery  the  collision  time  is  much  larger  than  the  pinch  time, 
and  the  distribution  function  is  no  longer  in  equilibrium.  Instead,  the  distribution 
function  changes  on  a  timescale  much  faster  than  ree.  Inelastic  processes  can  be 
ignored,  but  the  magnetic  field  is  sufficiently  strong  to  avoid  electron  runaway,  i.e. 
P  3>  1.  The  opposite  case,  electron  runaway,  is  an  old  but  still  active  research 
topic  that  is  not  at  issue  here. 

A  numerical  solution  with  parameters  P2  =  l  and  E2  !E\  —  \  was  done  to  see  the 
effect  of  heating.  Initially  the  distribution  was  cold  enough  for  the  collision  time 
to  be  comparable  or  smaller  than  the  timescale  for  heating,  and  the  distribution 
remained  close  to  maxwellian.  Eventually  the  plasma  heated  up  enough  for  the 
collision  time  to  become  large  compared  to  the  heating  time.  The  temperature 
increases  at  an  ever  decreasing  rate. 

Without  collisions  the  heating  can  be  calculated  analytically.  In  this  case  the 
relevant  equation  for  f()  becomes,  setting  C  and  D  to  zero. 

Oh  =  (E/P.E0)2  0  {  e  Oh  | 

Or  £2  0(, 1 P~2  +  Ov 

A  self-similar  solution  is  possible  for  P£3  1.  i.e.  excluding  the  low-velocity 

regime  where  the  abridged  equation  is  not  valid  anyway.  Ignoring  P-2  in  the  de¬ 
nominator  the  self-similar  solution  is28'29 

h  -x  so  *  exp(— r  /Yfo).  (45 1 

where  £0P)  increases  with  time  according  to 

=  £o(0)  +  5  f  )*  dt. 

JQ  toll 
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(46) 


When  the  exp(— £5  )-dependence  is  extended  to  small  £,  conservation  of  elec¬ 
trons  gives  the  factor  £q3.  The  energy  per  particle,  /  £4h  d£  =  £2  =  T.  the 
temperature,  is  then  increasing  in  time  more  slowly  than  linear.  £'2  x  t2^5 .  The 
(normalized)  power  input  per  particle  decreases  with  time  according  to 


9r3/s  oc  r-3/2, 


(47 


the  opposite  temperature  dependence  than  found  for  the  Spitzer  conductivity.  The 
conductivity  as  defined  here  is  <r=(power  input/E2),  which  is  3£2  jdt.  since  the  heat 
capacity  is  constant. 

That  the  temperature  dependence  of  the  conductivity  differs  from  the  usual 
result  should  not  surprise:  the  dominant  effects  differ.  The  Spitzer  conductivity  is 
found  from  a  balance  between  the  acceleration  by  the  electric  field  and  the  collision 
time,  assuming  the  temperature  remains  constant,  while  in  contrast  here  the  accel¬ 
eration  in  a  weakly  collisional  but  highly  magnetized  plasma  is  balanced  with  the 
temporal  development  of  the  (non-equilibrium)  distribution  function. 

The  electrical  conductivity  in  the  plasma  periphery  used  in  the  hydrodynamics 
computations  should  be  replaced  by  the  formulas  obtained  here  as  the  changes  in 
density  and  temperature  in  the  periphery  have  brought  the  plasma  in  the  weakly 
collisional  regime.  A  reasonable  transition  point  is  when  the  collison  time  is  1-10 
x  the  macroscopic  timescale  in  the  local  plasma,  e.g  dt/dinT  for  the  tempera¬ 
ture  increase.  When  the  transition  temperature  Tc  is  reached  in  the  MHD  codes 
the  available  conductivity  can  easily  be  modified  by  simply  replacing  the  tempera¬ 
ture  dependence  from  (T/Tc)3^2  to  (T/Tc)-3^2 ,  while  keeping  the  constant  at  the 
transition  temperature  unchanged.  Alternatively,  the  expression  derived  here  can 
be  used  once  it  is  smaller  than  the  Braginskii  conductivity.  However,  note  that 
( E/E0fl )2  =  (v0 /vq)2  x  Zj  1  2tt  <C  1,  with  vD  =E/B  the  drift  velocity. 

In  the  corona  of  laser-generated  plasma  measurements  of  the  electron  distribu¬ 
tion  function30  show  a  flat-topped  distribution  in  the  energy  regime  accessible  by 
the  diagnostics,  in  agreement  with  the  theory26. 


6.4  INELASTIC  COLLISIONS. 

This  section  estimates  the  decrease  in  ionization  rate  from  the  reduction  in  the 
distribution  function's  tail,  with  parameters  relevant  to  the  pinch  periphery.  From 
Table  1  it  is  apparent  that  a  stationary  distribution  should  exist  on  a  nano-second 
timescale.  The  stationary  distribution  thus  satisfies  the  equation 


Here  the  coefficients  7(e)  =  D (£)=/'  tl/2  f  dt. 
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When  the  normalized  electron  energy  t  is  large  compared  to  1  the  coefficient 
7.  the  (incomplete)  number  density,  goes  to  1:  the  coefficient  3.  the  (incomplete) 
energy,  also  asymptotes  to  1.  For  a  maxwellian  with  average  energy  3e0/2.  /  -x 
exp(— e/e0),  where  7e0  =  3. 

The  distribution  function  has  to  satisfy  a  number  of  conditions.  At  the  bound¬ 
ary  e  — +  00  f  and  df  /dt  are  zero,  but  there  is  no  such  condition  at  e  =  0.  Instead, 
there  is  the  important  condition  that  f  remain  positive,  which  is  trivially  satisfied, 
and  the  normalization  condition,  /0°°  A2/  dt  =  l.A  second  condition  may  follow 
from  physical  requirements,  e.g.  one  may  wish  to  have  the  energy  loss  in  the  inelas¬ 
tic  collisions  compensated  exactly  by  an  explicit  energy  source,  such  as  the  electric 
field. 


A  numerical  solution  is  straightforward  provided  the  equations  can  be  lin¬ 
earized,  i.e.  a  previously  calculated  distribution  function  is  used  to  find  the  co¬ 
efficients  7  and  3-  This  is  a  reasonable  procedure  because  the  coefficients  are  in¬ 
sensitive  to  the  exact  shape  of  the  distribution.  In  doing  this  another  requirement 
may  be  that  the  temperature  of  the  plasma  correspond  to  the  temperature  of  the 
distribution  that  enters  the  coefficients.  Finally,  the  solution  can  be  iterated  until 
convergence  is  obtained. 

An  interesting  feature  of  the  solution  are  the  large  changes  in  the  tail,  while 
the  body  remains  the  same.  Now  the  question  is,  what  do  we  want  with  this? 

One  especially  interesting  case  occurs  when  the  inelastic  collision  rate  is  large 
(compared  to  the  electron-electron  scattering  rate)  for  a  process  at  moderate  low 
energy  that  in  itself  is  not  critical.  The  distribution  function's  tail  will  then  be 
substantially  affected,  and  the  collision  rate  is  suppressed  at  a  higher  energy  where 
the  interesting  process  takes  place.  Something  similar  should  happen  when  all 
the  inelastic  processes  are  taken  into  account.  Estimates  along  the  same  lines  as 
presented  here  are  possible,  but  will  not  be  highlighted  because  they  would  be 
more  believable  with  a  numerical  code,  whose  development  is  already  sufficiently 
motivated. 


6.5  ELECTRIC  FIELD  EFFECTS. 


An  equilibrium  distribution  in  the  absence  of  inelastic  collisions,  but  with  an 
electric  and  magnetic  field  is  formally 


exp 
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with  the  additional  condition  that  the  normalisation  be  maintained. 


/ 


d£  ?h  =  1. 


(51) 


For  a  weak  electric  field  E  the  asymptotic  forms  of  C  and  3£D  are  unity,  and 
the  distribution  reduces  to  the  equilibrium  Maxwellian  h  ~  exp(—  3£2/2).  When 
the  magnetic  field  is  large,  such  that  ft£3  <C  1  over  most  of  the  interesting  region, 
the  electric  field  term  becomes  a  constant, 


3^(  ft-2  +  £6 


In  this  limit  the  term  is  small  compared  to  unity,  basically  because  anywhere  away 
from  the  central  region  in  the  pinch  (uo/uo)2  <C  1.  Hence  the  electric  field  does 
not  modify  the  denominator  3 £D  ~  1,  and  the  electron  distribution  is  basically 
thermal. 


In  contrast,  along  the  pinch  axis  ft  =  0,  and  the  integrand  in  the  exponent 
converges,  showing  that  the  present  formulation  is  not  rigorous  in  a  runaway  regime. 

Even  though  the  integral  does  not  converge  the  expression  can  be  a  good  ap¬ 
proximation  to  the  distribution  function  provided  the  integral  is  large  compared 
to  1,  i.e.  the  exponential  factor  is  close  to  zero.  Then  the  term  ( E/E0 )2  in  the 
denominator  suggests  an  enhancement  to  the  electron  temperature. 

In  any  case,  it  seems  that  effectively  hot  electron  distributions  can  be  generated 
with  electric  fields,  but  only  along  the  pinch  axis  where  the  magnetic  field  vanishes. 
Let  the  boundary  for  this  effect  be  uiceTe,  =  1  (ft  =  1).  In  numerical  units,  this 
occurs  at  a  radius  r  inside  the  pinch  where 


r(mm) 


_ n(cm3) _ 

3  x  1018  I (XI A)  Z  T3'2(eV )  ' 


(53) 


assuming  that  the  current  density  is  constant  over  the  pinch.  Thus,  this  argument 
does  not  apply  to  the  runin  stage,  but  only  to  the  stagnation  stage  of  the  pinch 
already  collapsed  on  axis. 

We  must  now  ask  if  there  is  experimental  evidence  for  tighter  bright  spots 
when  Z  is  large.  Riordan31  compares  the  emission  from  pinches  with  different  ma¬ 
terials.  These  indeed  show  that  tungsten  has  smaller  bright  spots  than  aluminum, 
intriguing.  Further  work  might  shed  clearer  light  on  such  phenomena. 
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TABLE  I.  Symbol,  edge  and  periphery. 


Symbol 

edge 

periphery 

Current 

I 

2  MA 

Radius 

r 

1  mm 

20  mm 

Magnetic  field 

B 

1000  T 

50  T 

Gyrofrequency 

oT 

II 

3 

1014/s 

5xl012 

gyroradius 

re 

10-5cm 

10~5  cm(  ? ) 

plasma  frequency 

“p 

1015/s 

10t2/s(?) 

Electron  density 

ne 

lO20 /cm3 

1014/cm3(?) 

Ion  density 

Tlt 

1019/cm3 

10 13 /cm3  (?) 

Typical  ionic  charge 

Z 

10 

10 

ionic  size 

a0/Z 

10-9cm 

Electron  temperature 

Te 

300  eV 

10  keV  (?) 

Ion  temperature 

T, 

300  eV  (?) 

1  keV  (?) 

Debye  length 

10-6cm 

3xl0-5cm 

Plasma  parameter 

ne  Aj, 

100 

3 

e-e  collision  time 

Tec 

1.5  ps 

100//s(?) 

e-i 

Tei  =T„/Z 

0.15  ps 

10/xs(?) 

e-i  energy  exchange 

Tt  =  ( M/m)Tet 

30  ns 

ms  (?) 

collisionality 

Uce  Tei 

30  >  1 

107  »  1(?) 

ionization:  Ar 

energy  2  keV 

rate  coefficient 

k  =<  av  > 

10-u  cm3/ s 

inelastic  collision  time 

rc 

10  ns 

Charge/length 

Ne 

6  Coul/cm 

drift  velocity 

VD 

3  x  lQ5cm/s 

drift  distance  in  rc 

3  x  lO"3 cm 
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